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CALCULATION  OF  COMPONENT  FORCES  AND  MOMENTS  OF  ARBITRARILY 
BANKED  CRUCIFORM  MISSILES  WITH  CONTROL  DEFLECTIONS 

by  Michael  J.  Hemsch,  Charles  A.  Smith, 

Jack  N.  Nielsen,  and  Stanley  C.  Perkins,  Jr. 

Nielsen  Engineering  & Research,  Inc, 

1.  INTRODUCTION 

Under  the  sponsorship  of  the  Office  of  Naval  Research  (ONR) , Nielsen 
Engineering  & Research,  Inc.  (NEAR)  has  been  developing  computer  programs 
for  the  prediction  of  the  aerodynamics  of  arbitrarily  banked  cruciform 
missiles  at  high  angles  of  attack.  The  first  two  years  of  the  work  were 
devoted  to  the  development  of  a computer  code  based  on  linear  supersonic 
flow  theory  for  the  calculation  of  loadings  on  cruciform  wing-body 
combinations.  The  results  are  reported  in  references  1 and  2.  In  1974- 
1975,  under  the  sponsorship  of  the  Naval  Weapons  Center  (NWC) , China 
Lake,  CA,  NEAR  developed  an  engineering  method  for  calculating  the  induced 
rolling  moments  of  cruciform  missiles  at  angles  of  attack  up  to  20°.  At 
the  same  time,  under  the  joint  sponsorship  of  NWC  and  NASA/Ames  Research 
Center,  NEAR  obtained  a systematic  set  of  component -buildup  data  on  a 
cruciform  body-cc.nard-tail  missile  up  to  24°  angle  of  attack  at  transonic 
and  supersonic  speeds.  The  data  included  individual  fin  loads.  The 
engineering  method  is  reported  in  reference  3 . The  data  are  presented  in 
references  4 and  5. 

Under  the  co-sponsorship  of  ONR  and  NASA/ARC,  NEAR  is  continuing 
and  extending  the  work  reported  in  reference  3.  The  work  involves  two 
tasks;  (1)  development  of  a single  comprehensive  code  for  the  engineering 
calculation  of  all  missile  component  aerodynamic  forces  and  uoments, 
except  minimum  drag,  for  angles  of  attack  up  to  the  unsteady  regime;  and 
(2)  comprehensive  comparison  of  predicted  results  with  systematic  data 
from  a variety  of  configurations.  The  work  is  to  be  accomplished  in  two 
phases.  The  first  phase  includes  development  of  a preliminary  version 
of  the  new  code  based  on  the  method  of  reference  3.  The  new  code  is  then 
to  be  compared  with  available  systematic  data  up  to  approximately  20° 
angle  of  attack.  This  report  describes  that  work.  The  work  planned  for 
the  second  phase  is;  (1)  correction  of  deficiencies  in  the  method  which 
were  uncovered  in  the  first  phase  through  comparison  of  experiment  and 
theory,  (2)  extension  of  the  method  to  angles  of  attack  up  to  the  unsteady 
regime,  and  (3)  comparison  of  the  extended  method  with  systematic  high 
angle  of  attack  data  obtained  for  that  purpose. 
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When  the  work  reported  herein  was  in  the  planning  stage,  it  was 
thought  that  in  the  first  phase  a comprehensive  comparison  between 
theoretical  predictions  and  data  for  a number  of  different  missile 
configurations  and  Mach  numbers  would  be  appropriate.  However,  detailed 
comparisons  of  theory  and  experiment  for  individual  fin  loads  were  found 
to  be  essential  for  uncovering  and  sorting  out  the  many  and  varied 
phenomena  involved.  As  a result,  the  comparisons  are  directed  principally 
toward  systematic  data  which  include  individual  fin  loads.  Where  discrep- 
ancies between  theory  and  data  are  found,  remedies  are  suggested  and 
further  work  outlined.  To  insure  that  the  phenomena  uncovered  are  not 
the  result  of  isolated  cases,  several  different  configurations  are 
studied. 

The  engineering  method  to  be  described  in  this  report  is  composed 
of  five  basic  procedures:  (1)  computation  of  forces  and  moments,  with 

no  external  vortices  present,  for  individual  fins  in  the  nresence  of  a 
body  at  arbitrary  angle  of  attack  and  bank  angle  with  control  deflections; 
(2)  computation  of  forces  and  moments  induced  on  individual  fins  due 
to  the  presence  of  external  vortices;  (3)  calculation,  using  an  empirical 
method,  of  the  strengths  and  positions  of  body  vortices  generated  upstream 
of  the  first  set  of  fins;  (4)  computation,  using  a vortex  tracking  method, 
of  the  strengths  and  positions  of  vertices  between  the  first  and  second 
sets  of  fins;  and  (5)  determination  of  body  forces  and  moments  due  to 
the  presence  of  the  fins  and  vortices. 

This  report  and  reference  3 are  intended  to  be  sufficient  for  an 
understanding  of  the  method.  Reference  3 gives  a detailed  description 
of  procedures  (1)  and  (3) . Although  detailed  descriptions  of  procedures 
(2)  and  (4)  are  available  in  reports  published  prior  to  reference  3,  for 
the  sake  of  completeness,  that  work  is  made  available  in  this  report 
in  Appendices  A through  D.  Section  2 of  this  report  presents  an  overall 
introduction  to  the  method  and  describes  procedure  (5),  the  development 
of  which  is  Task  1 of  the  present  effort.  Detailed  comparisons  between 
data  and  theory  for  determining  the  effectiveness  of  the  above  procedures 
(Task  2)  are  presented  in  sections  3 to  5.  Section  3 is  concerned  with 
cases  in  which  external  vortex  effects  are  unimportant  (procedure  (1) ) . 
Section  4 considers  situations  for  which  the  effects  of  body  vortices 
are  important  (procedures  (2)  and  (3))  but  do  not  require  the  vortex 
tracking  of  procedure  (4) . Section  5 compares  theory  and  data  for  cases 
in  which  all  of  the  above  procedures  are  required.  A user's  manual  for 
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the  Task  1 code  is  presented  in  Appendix  E.  Appendix  F contains  exten- 
sions made  to  a computer  code  developed  earlier  by  NEAR  for  ONR  (refs.  1 
2,  and  3) . The  additions  were  necessary  for  the  computation  of  certain 
interference  quantities  used  for  the  calculation  of  body  forces  in  the 
presence  of  cruciform  fins.  Appendix  F also  includes  a study  of  optimuni 
panel  spacing  to  be  used  for  that  program. 


2.  INTRODUCTION  TO  THE  ENGINEERING  METHOD 
AND  DESCRIPTION  OF  THE  EXTENSIONS  FOR 
CALCULATING  BODY  FORCES  AND  MOMENTS 


2.1  Introduction  to  the 
Engineering  Method 

In  developing  the  computer  program,  we  have  adopted  an  analytical 
approach  to  the  maximum  extent  possible  compatible  with  obtaining  a 
complete  program  for  all  forces  and  moments  acting  on  the  missile  com- 
ponents except  for  minimum  drag.  However,  it  has  been  necessary  to 
adopt  semi-empirical  techniques  for  parts  of  the  method.  By  comparisons 
of  systematic  data  with  theoretical  predictions,  the  shortcomings  and 
limitations  can  be  determined  and  necessary  improvements  identified. 

To  place  the  present  method  in  perspective,  a brief  description  of 
the  phenomena  involved  in  flow  over  cruciform  missiles  is  presented. 
Consider  the  body -canard-tail  configuration  shown  in  figure  1.  Experi- 
mental work  on  such  configurations  has  shown  that  the  following  sources 
of  component  interference  can  be  important  for  included  angles  of  attack 
in  the  range  of  0 < a ^ 20°. 


1.  Panel-panel  interference;  Consider  a set  of  opposing  fins  at 
angle  of  attack  and  roll  joined  together  at  their  root  chords.  We  will 
call  such  a configuration  a "wing  alone."  A body  can  be  inserted  between 
the  fins  and  another  set  of  fins  can  be  added  perpendicular  to  the  first 
set.  Any  change  in  loading  experienced  by  a fin  due  to  changes  in  loading 
of  another  fin  is  said  to  be  due  to  panel-panel  interference  and  is 
affected  by  the  presence  of  the  body. 

2.  Nose  vortex-fin  interference:  At  sufficiently  high  angles  of 

attack,  separated  flow  will  develop  on  the  missile  upstream  of  the  first 
set  of  fins.  The  presence  of  the  separated  flow  over  the  fins  will  alter 
their  loading. 

3.  Canard  vortex-tail  interference:  Each  of  the  fins  in  the  first 

finned  section  will  shed  trailing  vorticity.  That  vorticity  passing  the 
tail  section  will  alter  the  loading  experienced  by  the  tail  fins. 

4.  Afterbody  vortex-tail  interference:  At  sufficiently  high  angles 

of  attack,  the  afterbody  boundary  layer  will  separate.  The  shed  vorticity 
will  roll  up  and  interact  with  the  tail  fins  altering  the  loadings  they 
will  experience. 
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As  discussed  in  the  introduction,  the  present  method  models  the 
above  phenomena  with  five  basic  procedures.  The  first  procedure  computes 
fin  forces  and  moments  when  no  external  vortices  are  present;  that  is, 
no  vorticity  generated  upstream  of  the  fins  is  present.  To  do  this  we 
have  developed  a correlation  method  based  on  the  concept  of  an  "equivalent 
angle  of  attack."  The  equivalent  angle  of  attack,  agq»  is  defined  to  be 
the  angle  of  attack  of  a wing  alone  such  that  the  normal  force  of  the 
wing  alone  is  twice*  the  actual  normal  force  on  a fin  mounted  on  a body 
with  arbitrary  deflection  and  roll  orientation.  The  various  interference 
factors  which  account  for  fin  deflection,  roll  orientation  and  ratio  of 
fin  span  to  body  radius  are  described  in  reference  3.  The  Task  1 computer 
code  requires  that  the  wing-alone  characteristics  of  normal  force  and 
center  of  pressure  location  be  input  as  functions  of  angle  of  attack.  The 
usefulness  of  the  a^q  concept  is  demonstrated  in  reference  6,  pages  21 
to  25.  The  present  version  of  this  procedure  is  limited  to  symmetric 
deflection  of  opposing  fins.  Furthermore,  only  the  first  set  of  fins 
may  be  deflected.  Both  of  these  limitations  are  easily  removed  and  may 
be  during  the  second  phase  of  the  work.  The  fins  must  not  be  cambered 
or  twisted. 

The  second  procedure  mentioned  above  is  concerned  with  the  changes 
in  fin  forces  and  moments  when  vortices  generated  upstream  of  the  fins 
are  present.  A full  description  of  the  method  used  is  given  in 
Appendices  A and  B and  will  not  be  repeated  here.  However,  the  reader 
should  be  aware  that  the  procedure  assumes  that  the  flow  induced  by  the 
vortices  can  be  described  by  linear  theory  and  uses  reverse  flow  methods 
to  obtain  the  appropriate  equivalent  angle  of  attack  increment  due  to 
the  vortices.  Because  of  the  flow  modeling  used  to  track  upstream 
vortices  over  a finned  section,  some  restrictions  on  fin  leading-edge  and 
trailing-edge  sweep  exist.  See  Appendix  A for  further  discussion. 

For  the  third  procedure  the  separated  flow  over  the  body  upstream 
of  the  first  set  of  firs  is  modeled  as  two  symmetric  potential  line 
vortices.  The  body  is  assumed  to  be  axisymmetric.  It  is  also  assumed 
to  be  cylindrical  downstream  of  the  nose.  The  missile  station  at  which 
separation  begins,  thac  is,  the  station  at  which  the  vorticity  is  rolled 
up,  is  estimated  from  an  empirical  correlation.  The  vortex  positions  on 
the  nose  are  also  obtained  from  an  empirical  correlation.  The  strengths 


*Recall  that  the  wing  alone  is  composed  of  two  opposing  fins. 
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of  the  nose  vortices  at  the  leading  edge  of  the  first  finned  section  are 
then  calculated  using  the  vortex  impulse  theorem  and  two-dimensional 
crossflow  drag  coefficient  correlations.  A full  description  of  this 
procedure  is  given  in  reference  3. 

The  fourth  procedure  is  concerned  with  the  wake  over  the  afterbody 
section.  The  afterbody  is  assumed  to  be  cylindrical.  In  this  region, 
all  the  vortices  present  are  tracked  using  slender-body  theory  from  the 
trailing  edge  of  the  first  set  of  fins  to  the  leading  edge  of  the  second 
set  of  fins.  Appendix  C describes  the  method  used  to  obtain  the  initial 
strengths  and  positions  of  the  vortices  shed  by  the  first  set  of  fins. 
Criteria  for  determining  the  axial  location  of  the  onset  of  boundary- 
layer  separation  and  roll  up  of  the  shed  vorticity  from  the  afterbody 
are  based  on  the  local  crossflow  velocities  and  the  presence  or  absence 
of  nose  vortices.  If  afterbody  vortices  form,  they  are  assumed  to  grow 
as  two  potential  line  vortices  under  the  influence  of  the  vorticity  shed 
upstreaun.  They  are  allowed  to  grow  unsymmetrically  if  the  local  flow  is 
unsymmetrical  about  the  plane  containing  the  body  and  wind  aixes.  The 
strengths  of  the  afterbody  vortices  are  assumed  to  be  given  by  the  vortex 
impulse  theorem  and  two-dimensional  crossflow  drag  correlations.  A 
complete  description  of  this  procedure  is  given  in  Appendix  D. 

Procedure  (5)  which  is  concerned  with  body  forces  and  moments  is 
described  in  the  remainder  of  section  2. 


2.2  Description  of  Extensions 
to  the  Method 


2.2.1  Nose  section  forces  and  moments.-  The  present  version  of  the 
method  assumes  that  the  nose  vortices  are  symmetrical.  Hence,  in  unrolled 
body  coordinates,  the  nose  side  force,  yawing  moment  and  rolling  moment 
are  zero  (see  fig.  2 for  a description  of  the  coordinate  system  used  in 
this  report) . The  normal-force  coefficient  is  determined  from  the 
expression 


'M  “ 


+ c sin^a^  ^ 
C R 


(1) 


If  the  user  does  not  input  Cjj  , the  program  uses  the  slender-body  theory 
result,  Sa  “ 2.  The  two-dimensional  crossflow  drag  coefficient  is 
obtained  from  Fidler's  correlations  (ref.  7)  for  crossflow  Mach  numbers, 
M^,  less  than  0.7  and  from  Jorgensen's  correlation  (ref.  8)  for 
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(a)  Total  loads  and  fin  numbering  system. 
Figure  2.  Definition  of  coordinate  systems 


Yaw  plane 


(b)  Sign  convention  for  canard  deflection  angles 


Figure  2.  Continued. 


CBM2, 

CRM2 


looking 

forward 


CBM4, 

CRM4 


JOTE 

(c) 


Hinge 

line 


Normal  forces  are  measured  perpendicular  to  the 
fin  planform  and  are  considered  to  act  at  the 
intersection  of  the  fin  hinge  line  and  root  chord. 


Axis  systems  and  positive  sign  convention  for  fins. 
Figure  2.  Concluded. 
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Mj,  > 0.7.  Note  that  the  usual  formulation  of  equation  (1)  includes  a 
multiplicative  factor  for  finite  length  effects  which  is  not  used  in  the 
present  method. 

The  center  of  pressure  position  for  the  linear  term  of  equation  (1) 
is  obtained  from  the  slender-body  theory  expression 


^linear 


(2) 


where  Xjj  is  the  length  of  the  nose  and  Vjj  is  the  nose  volume.  The 
code  obtains  and  Vn  by  numerically  integrating  the  expressions: 


S 


C 


2 


r (x)  dx 


(3) 


V, 


N 


r^  (x)dx 


(4) 


using  the  input  nose  coordinates.  The  quantity  x^  is  the  estimated 
initial  axial  position  of  flow  separation  on  the  nose  (see  ref.  3, 
pp.  10-11) . Since  the  crossflow  drag  coefficient  as  used  in  the  present 
method  is  not  a function  of  axial  position  on  the  body,  the  axial  center 
of  pressure  location  of  the  viscous  portion  of  the  nose  normal  force 
(second  term  of  eq.  (1))  is  given  by 


X 


viscous 


r (x)  dx 


(5) 


In  reference  3,  two  methods  for  computing  nose  vortex  strengths 
were  used  depending  upon  the  length  of  the  body  up  to  the  leading  edge 
of  the  first  finned  section.  During  the  course  of  the  work  for  Phase  I 
(reported  herein) , it  was  found  that  the  method  based  on  the  crossflow 
drag  coefficient  was  sufficient.  The  code  uses  only  that  method  (see 
ref.  3,  pp.  12-13). 

2.2.2  Finned  section  forces  and  moments.-  In  reference  3,  a method 
was  presented  for  computing  the  normal-force  coefficient  for  a rolled 
fin  at  angle  of  attack  including  the  effects  of  panel-panel  interference 
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and  vortex-fin  interference  (see  ref.  3,  pp.  14-24  and  Appendices  A 
and  P of  this  report) . A method  for  computing  the  fin  rolling -moment 
coefficient  was  also  given  (see  ref.  3,  pp.  29-30) . This  section 
describes  the  computation  of  all  fin  forces  and  moments. 

Once  the  equivalent  angle  of  attack  for  a given  fin  has  been 
computed,  the  fin  normal-force  coefficient  and  center  of  pressure  loca- 
tion can  be  found.  The  rolling  moment  for  fin  i has  been  given  in 
reference  3 as 

CRMi  = cos  |cNi(agq^p)  ■ ^ 

* [°*K%q,P  * ““eq>v)  ' 

where  is  the  equivalent  angle  of  attack  without  the  vortices 

presc'.*  and  (ACg^)^  is  the  angle  of  attack  induced  on  the  fin  by  the 
vor^'f'.  s.  The  fin  bending  moment  at  the  root  chord  follows  from  the 
normal  force  and  rolling  moment 

CBMi  = ■ ■ - CNi  (7 

cos  0 ■ £ 

1 r 

The  hinge  moment  is  given  by 

^ 

The  contributions  of  the  fin  forces  and  moments  to  the  total 
missile  forces  and  moments  are  given  next  in  unrolled  body  coordinates. 

C - (CN2  cos  + CN4  cos  6^)  cos  (J 

“fins 

+ (CNl  cos  6,  + CN3  6_) sin  0 (9 

1 


CL.  “ (CN2  cos  6_  + CN4  cos  6 ) sin  (f> 
^fins  ^ 

- (CNl  COS  6^  + CN3  cos  6^)  cos  0 


(11) 


to. 


2.  » CNl  sin  6,  + CN2  sin  6 + CN3  sin  6 + CN4  sin  6 

*finS  12  3 4 


m 


fins 


" NfinsV  ) 


+ (CBM2  sin  6g  - CBM4  sin  6^  + CHMl  + CHM3) sin  0 

+ (CBMl  sin  6,  - CBM3  sin  6„  + CHM2  + CHM4) cos  0 
1 


n 


fins 


■fins 


(12) 


- (CBM2  sin  6 - CBM4  sin  6,  + CHMl  + CHM3)cos  0 


+ (CBMl  sin  6^  - CBM3  sin  63  + CHM2  + CHM4) sin  0 


C..  = - CRMl  + CRM2  + CRM3  - CRM4 

^fins 


The  total  lift  and  drag  due  to  lift  for  the  fins  are 


'T  “ sin  a 

^fins  Nfins  ^fins  ^ 


C = C sin  a + C.  cos  a 

°fins  Nfinf  ^ Afins 


(13) 

(M) 

(15) 

(16) 


In  reference  9,  a method  based  on  slender-body  theory  is  presented 
for  computing  the  interference  lift  produced  on  a body  at  zero  bank  in 
the  presence  of  a set  of  all-movable  fins.  For  the  case  of  ^ 0 and 
6 = 0,  an  interference  factor,  Kg,  is  defined  such  that 


dc, 


’N, 


'N, 


= K 


W 


B (W) 


B da 


* a. 


6=0 


a=0 


(17) 


The  normal-force  coefficient  produced  on  the  fins  in  the  presence  of  the 
body  is  given  by 


dC. 


N, 


w 


a“0 


• a , 

c ’ 


6 = 0 


(18) 


il 
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Dividing  (17)  by  (18)  gives 


'N, 


B(W) 


"N, 


W(B) 


(19) 


For  the  case  of 
such  that 


a = 0 and  6^0.  an  interference  factor,  k_,  is  defined 

C B 


<3C, 


N, 


W 


'N, 


B (W) 


B da 


a_  = 0 


a»0 


(20) 


The  normal-force  coefficient  produced  on  the  fins  in  the  presence  of  the 
body  is  given  by 


= k W 

'Nw(B)  ^ 


. 6 


cp*0 


a 


(21) 


Dividing  (20)  by  (21)  gives 


”b(W)  _ 


(22) 


Nielsen  (ref.  9)  has  shown  that  the  ratios  and  differ  by 

less  than  5 percent  over  the  entire  range  of  the  radius  to  semispan  ratio. 
This  result  indicates  that  the  fraction  of  wing  normal  force  which  is 
carried  over  onto  the  body  can  be  assumed  to  be  independent  of  the  means 
of  lift  production.  Hence,  the  present  method  assumes 


Similarly, 


c 

Nbody  .^3 


Sr  “ ^ ^ 

*body  *fins 


(25) 


(24) 


Reference  10  gives  a method  for  computing  the  axial  location  of 

the  center  of  pressure,  x , for  the  force  on  the  body  due  to  the  fins 

B 
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for  zero  fin  deflection  and  zero  roll.  We  assume  for  the  present  method 
that  Xg  gives  the  correct  value  independent  of  the  means  of  lift 
production  on  the  fins.  The  moments  created  by  the  body  force  due  to 
the  presence  of  the  fins  are 


C = C 

"’body 


body 


'4. 


^MC  ~ 


’^odv 


body 


C.  =0  (27) 

■*body 

The  user  has  the  option  of  having  the  code  compute  Kg  and  Xg 
using  the  methods  of  reference  10  or  of  using  values  obtained  by  some 
other  method  such  as  the  code  of  Appendix  F. 

2.2.3  Afterbody  section  forces  and  moments.-  The  model  used  in  the 
present  method  for  vortex  tracking  over  the  afterbody  section  has  been 
discussed  previously  in  reference  3,  pages  24-28,  and  in  reference  11, 
pages  6-11.  The  discussion  is  repeated  and  extended  in  Appendix  D.  The 
present  code  does  not  account  for  tearing  of  afterbody  vortices. 

The  afterbody  loading  is  computed  from  the  vortex  impulse  theorem 
using  the  vortex  strengths  and  positions  computed  by  the  -cracking  method. 
The  vortex  impulse  theorem  for  moderate  angles  of  attack  gives  the 
following  expression  for  the  loading  on  an  infinitesmal  length  of  the 
afterbody  in  the  presence  of  vortices 

NV(x) 

‘ ^ E ■’‘'■'a')  (28) 


where  a*  is  the  nondimens ional  complex  distance  between  a vortex  and 
its  image  inside  the  cylinder.  Since  the  afterbody  vortices  may  form 
anywhere  between  the  canard  trailing  edge  and  the  tail  leading  edge, 
the  number  of  vortices  NV  is  a function  of  cixial  position. 


We  can  integrate  equation  (28)  to  obtain  the  afterbody  forces  and 
moments.  The  results  are 


- iC, 


4Tra‘= 
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XCTE. 


(29) 


NV(x)  XTLE 


C - iC  = 


4Tra‘" 
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£ I 


(Xj^C  - x)d(r'a') 


(30) 


j=l  XCTE 


The  integral  in  equation  (30)  is  approximated  by  numerical  quadrature  in 
the  code.  The  lower  and  upper  limits  of  the  integral  are  the  trailing 
edge  of  the  canard  section  and  the  leading  edge  of  the  tail  section 
respect ively . 
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3.  COMPARISONS  WITH  DATA  FOR  FINS  WITH 
NO  EXTERNAL  VORTICES  PRESENT 

Ilie  purpose  of  this  section  is  to  evaluate  the  ability  of  the  method 
to  predict  fin  loads  in  the  absence  of  external  vortices.  The  fins  and 
fin-body  combinations  to  be  considered  are  shown  in  figure  3.  Section  3.1 
considers  clipped  delta  fins  with  a leading-edge  sweep  of  45°  and  a 
taper  ratio  of  0.06  (aspect  ratio  = 3.3).  According  to  reference  12, 
such  fins  should  have  very  little  leading-edge  vortex  lift.  Hence, 
section  3.1  is  an  evaluation  of  the  method  for  fins  with  no  significant 
vortex  lift.  Section  3.2  considers  two  sets  of  low-aspect-ratio  fins 
for  which  fin-generated  vortex  lift  is  a substantial  fraction  of  the 
total  fin  lift. 

3.1  Moderate-Aspect-Ratio  Fins 

Data  for  this  set  of  fins  (C^)  mounted  as  shown  in  figure  3 (a)  were 
obtained  for  several  Mach  numbers,  bank  angles,  and  fin  deflections  for 
included  angles  of  attack  up  to  24°  (ref.  4) . The  tests  were  conducted 
in  such  a way  that  the  effects  of  bank  only,  deflection  only,  and 
combinations  of  both  could  be  studied.  A small  eimount  of  nose  vorticity 
was  observed  to  be  present  (ref.  5)  that  may  slightly  affect  the  fin  1 
loads. 

The  wing-alone  characteristics  used  for  the  predictions  are  shown 
in  figure  4.  The  normal  force  and  lateral  location  of  center  of  pressure 
data  were  obtained  for  < 30°  from  reference  4.  The  wing-alone 
values  were  computed  from  fin  loads  measured  for  <t>  - 0,  no  yaw  control 
and  various  values  of  pitch  control.  A description  and  experimental 
verification  of  this  procedure  is  given  in  reference  6,  pages  21-25.  The 
results  were  extrapolated  to  higher  angles  of  attack  using  a procedure 
given  by  Aiello  (ref.  13)  . The  axial  location  of  the  fin  center  of 
pressure  was  estimated  as  follows.  Axial  locations  of  the  center  of 
pressure  for  several  sets  of  closely  related  fins  tested  alone  and  on 
bodies  (refs.  4,  14,  and  15)  were  compared  for  equivalent  angles  of  attack 
up  to  45°  (see  fig.  4) . A linear  best  fit  was  applied  to  the  data  wher- 
ever the  scatter  was  not  large.  For  M^  = 0.8,  the  scatter  was  too  large 
to  fit  for  < 14°  and  x/c^  was  arbitrarily  set  equal  to  0.55.  In 
each  of  the  sut»sections  below  the  discussion  centers  first  on  the  effects 
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NOTE;  All  dimensions  are  in  inches. 


(a)  Body-canard  configuration  no.  1 (BC  1) ; moderate-aspect- 
ratio  Cg  fins  mounted  on  3-caiiber  tangent  ogive  nose 
with  7-caliber  cylindrical  afterbody. 

(MICOM-NWC  tests,  ref.  4). 


(b)  Body-canard  configuration  no.  2 (BC  2);  low-aspect- 
ratio  fins  mounted  on  3-caliber  tangent  ogive 
nose  with  7-caliber  cylindrical  afterbody 
(MICCM-NWC  tests,  ref.  4). 


(c)  Body-tail  configuration  no.  1 (BT  1);  low-aspect- 
ratio  T fins  mounted  on  10-caliber  body 
(MICCM  tests,  ref.  16). 


Figure  3.  Configurations  used  for  comparisons 
of  data  and  theory. 
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1^ ' 
1.5  2.77 


(g)  Body-canard  configuration  no.  3 (BC  3) ; moderate- 
aspect-ratio  AIM-9L  canard  fins  mounted  on  20- 
caliber  cylindrical  afterbody  with  2-caliber 
nose  (AIM-9L  tests,  ref.  14)  . 


• V 
'i: 


(h)  Body-canard- tail  configuration  no.  1 (BCT  1) ; 
C and  T fins  moxinted  on  10-caliber  body. 
(MICOM-NWC  tests,  ref.  4). 


0 1.5  2.77  12.87  16.70 


(i)  Body-ceuiard-tail  configuration  no.  2 (BCT  2); 
AIM-9^j  model,  reference  14. 

Figure  3.  Concluded. 
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Fin  angle  of  attack,  a,  degrees 
(b)  Lateral  center  of  pressure  position 


Figure  4.  Continued 


of  bank  only,  second  on  the  effects  of  fin  deflection  only,  and  third  on 
the  combined  effects  of  bank  and  fin  deflection. 

3.1.1  Supersonic  flow.  Mco  *=  1.75.-  For  this  Mach  nximber  the  leading 
edges  of  the  fins  are  supersonic  for  = 0 and  0=0.  The  effects  of 
bank  only  are  presented  first  for  several  angles  of  attack.  Fin  normal 
forces  are  shown  in  figure  5 with  the  corresponding  fin  rolling  moments 
shown  in  figure  6.  The  agreement  is  good.  However,  a comparison  of  the 
data  for  fins  1 and  2 at  0 = 45°  reveals  some  discrepancy.  The  two 
fins  should  have  the  same  loads  at  0 = 45°.  Since  the  normal-force 
coefficients  of  fins  2 and  4 agree  at  0=0,  it  seems  likely  that  the 
data  for  fin  1 are  slightly  too  high. 

For  the  case  of  bank  with  no  fin  deflections  and  no  external  vortices, 
the  equivalent  angle  of  attack  for  each  fin  depends  only  on  the  and 

terms  of  equations  (23)  to  (26)  of  reference  3.  For  fins  1 and  2 at 
; = 45°,  the  fraction  of  the  total  fin  load  for  the  case  under  discussion 
represented  by  the  term  ranges  from  25  percent  at  = 12.5°  to 

40  percent  at  = 20.6°.  Thus,  the  generally  good  agreement  exhibited 
in  figure  5 suggests  that  the  present  method  adequately  accounts  for 
bank  angle  when  no  fin  vorticity  is  generated.  The  agreement  for  fin 
rolling  moment  is  similar  to  that  for  fin  normal  force. 

The  fxn  loads  with  fin  deflection  and  no  bank  are  shown  for 
“ 1.75  in  figure  7.  The  data  are  shown  as  open  symbols.  Predictions 
using  the  Task  1 code  are  shown  as  solid  and  dashed  lines.  The  agreement 
between  data  and  predictions  using  the  present  method  is  good  except 
for  fin  3.  For  fin  3 the  rate  of  change  of  normal  force  with  a^, 
although  small,  is  about  twice  as  great  as  predicted. 

The  combined  effects  of  fin  deflection  and  bank  are  shown  in  fig- 
ure 8.  The  agreement  is  about  what  one  would  expect  from  a comparison 
with  figures  5 and  7,  The  theory  assumes  that  fin  deflection  and  bank 
effects  cam  be  linearly  superposed  and  figures  5,  7,  and  8 seem  to  con- 
firm that  assumption.  Note  from  figures  5,  6,  and  7 that  the  behavior 
of  the  fin  rolling  moment  in  th=  absence  of  free  vortices  closely 
follows  that  of  the  fin  normal  force. 

3.1.2  Subsonic  Flow.  = O.8.-  At  high  angles  of  attack,  the  flow 
about  body-canard  configuration  no.  1 (see  fig.  3(a))  at  = 0.8  is 
transonic.  Hence,  the  results  of  this  section  represent  an  evaluation 
of  the  method  for  extreme  transonic  flow  conditions.  The  effects  of 


Figure  5 


0,  degrees 
(b)  Fin  2. 

Fin  normal  force  for  bank  with  no  fin  deflection 
moderate-aspect-ratio  fins  (Cg) , BC  1, 

- 1.75;  ■ model  base  area. 


34 


<p,  degrees 
(a)  Fin  1 


0,  degrees 
(b)  Fin  2. 

Figure  6.  Fin  rolling  moment  for  bank  with  no  fin 
deflection,  moderate-asp)ect-ratio  fins  (Cg)  , 

1,  M » 1.75;  Sp  - model  base  area, 

- model  base  diameter. 


Included  angle  of  attack,  a^,  degrees 
(a)  Normal-force  coefficient. 


Figure  7.  Fin  normal  force  and  rolling  moment  for  fin  deflection 
with  no  bank,  moderate- aspect-ratio  fins  (Cg) , 

BC  1,  = 1.75;  Sj^  = model  base  area, 

I,  * model  base  diameter. 


0 2 4 6 8 10  12  14  16  18  20  22  24 

Included  anfjle  of  attack,  a^,  degrees 

(b)  Rolling -moment  coefficient. 

Figure  7.  Concluded. 
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bank  only  are  presented  in  figure  9.  The  agreement  is  generally  very 
good  except  for  fin  1.  Again,  the  data  for  fin  1 and  fin  2 do  not  agree 
at  <p  ^ 0.  It  seems  likely  that,  for  this  case  too,  the  fin  1 data  are 
too  high.  Adjusting  the  fin  1 data  so  that  they  agree  with  the  fin  2 
data  at  0 “ 45°  gives  very  good  agreement  with  the  theory. 

The  fin  loads  due  to  fin  deflection  with  no  bank  are  shown  in  fig- 
ure 10.  The  agreement  is  good  except  for  fin  1.  Since  fin  1 should  have 
the  scune  load  at  ■=  0°  as  fin  3,  it  appears  that  the  data  for  at 
least  one  of  the  fins  are  in  error.  If  the  data  for  fin  1 are  adjusted 
so  as  to  agree  with  fin  3 at  = 0°,  it  is  found  that  the  theory 
slightly  overpredicts  the  decrease  in  fin  normal  force  as  a increases 
for  > 6°.  A possible  explanation  for  this  discrepancy  may  be  the 
change  in  fin  leading-edge  sweep  with  etc.  0=0  and  s o,  the 

leading  edge  of  fin  1 is  swept  45°.  As  increases,  the  effective 

sweep  of  the  leading  edge  increases.  It  seems  likely  that  the  fin  picks 
up  significant  vortex  lift  as  etc  increases.  This  point  will  be  dis- 
cussed further  in  section  3.2.  Note  that  for  > 16°,  fin  1 seems  to 
be  favorably  affected  by  the  nose  vorticity.  This  point  will  be 
discussed  further  in  section  4. 

The  combined  effects  of  fin  deflection  and  bank  are  shown  in  fig- 
ure 11.  Again,  the  agreement  is  good  except  for  fin  1.  The  discrepancy 
for  fin  1 is  about  what  would  be  expected  from  the  results  of  the  bank 
only  (fig.  9)  and  fin  deflection  only  (fig.  10)  comparisons  as  discussed 
above.  The  data  for  fin  deflection  only  are  suspect. 


3.2  Low-Aspect-Ratio  Fins 

3.2.1  Aspect  ratio  1.3  fins.-  Data  for  the  Cr  set  of  fins  (see 
fig.  3(b))  were  obtained  for  several  Mach  numbers  and  bank  angles  for 
model  angles  of  attack  up  to  24°  (ref.  4) . The  fins  were  not  deflected 
except  in  pitch  at  0=0°.  The  wing-alone  normal  force  was  computed 
using  the  same  procedure  used  for  the  moderate-aspect-ratio  fins  (see 
section  3.1) . The  cixial  and  lateral  locations  of  the  center  of  pressure 
were  estimated  by  the  methods  of  reference  13.  The  wing-alone 
characteristics  are  given  in  figure  12. 

Comparisons  of  theory  and  data  for  = 1.75  for  nose-canard 
configuration  no.  2 (C7  fins)  are  presented  in  figure  13.  The  data 
trends  are  similar  to  those  observed  in  figure  5 for  the  moderate-aspect- 
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Figure  11.  Fin  normal  force  for  combined  bank  and  fin 
deflection,  moderate-aspect-ratio  fins  (C  ) , BC  1, 
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Figure  11.  Concluded 


Fin  angle  of  attack,  a,  degrees 
(a)  Normal-force  coefficient. 


Figure  12.  Wing-alone  characteristics  for  low-aspect-ratio  fins  (C  and  T ) 
flR  ” 1-3,  S„  “ planform  area.  ''  ® 
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Fin  angle  of  attack,  a,  degrees 
(b)  Axial  center  of  pressure  position. 
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Fin  angle  of  attack,  a,  degrees 
(c)  Lateral  center  of  pressure  position 
Figure  12,  Concluded. 


ratio  fins  (Cg) . However,  the  predictions  are  low  for  fin  1 and  high  for 
fin  3.  Since  no  such  effect  was  observed  for  the  Cg  fins,  we  are  led 
to  attribute  the  anomaly  to  the  difference  in  fin  aspect  ratio.  The 
present  method  uses  the  small  angle  definitions  of  angle  of  attack  and 
angle  of  sideslip.  Hence,  for  fins  1 and  3,  we  have 

a = a^sin  0 (31) 

P = a^cos  0 (32) 


as  shown  in  the  sketch  below. 


Sketch  1.  Small  angle  definitions  of  fin  3 
angle  of  attack  and  angle  of  sideslip. 


(For  high  angle  of  attack,  these  definitions  must  be  modified.)  Note 
that 

I = cot  0 (33) 

y 

which  means  that  when  fins  1 and  3 are  undeflected,  the  fin  sideslip 
angle  is  greater  than  the  fin  angle  of  attack  for  0 < 45°.  For  the  low- 
aspect-ratio  fins  under  discussion  in  this  section,  it  seems  likely  that 
the  leading-edge  and  side-edge  vortex  lift  is  strongly  influenced  by  the 
sideslip  under  such  conditions.  The  theory  used  in  this  report  does  not 
take  this  phenomenon  into  account.  In  fact,  the  present  method  assumes 
that  the  ratio  of  vortex  lift  to  total  lift  at  a given  angle  of  attack 
does  not  change  with  angle  of  sideslip. 
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comparisons  of  theory  and  data  for  the  fins  for  = 0.8  are 

presented  in  figure  14.  Note  that  the  comparisons  for  fins  1 and  3 are 
similar  to  those  for  = 1.75.  The  theory  indicates  that  the  equivalent 
angles  of  attack  for  fins  1 and  3 are  less  than  20°.  Hence,  considering 
the  wing-alone  force  characteristics  of  figure  12,  we  can  assume  that  the 
flow  over  the  fins  is  essentially  attached  except  for  edge  vortices.  As 
discussed  above,  the  effect  of  sideslip  on  edge  vortex  lift  may  be  an 
important  phenomenon. 

For  discussion  purposes,  let  us  consider  a fin  with  aspect  ratio  1.3 
but  with  a delta  planform.  For  such  a fin,  all  the  vortex  lift  is  due  to 
the  leading  edge.  It  has  been  shown  for  such  fins  in  incompressible  flow 
at  a given  angle  of  attack  with  no  sideslip  that  the  ratio  of  vortex  lift 
to  attached  flow  lift  decreases  as  the  wing  leading-edge  sweep  decreases 
(ref.  12) . It  seems  likely  that  a similar  result  occurs  for  changes  in 
leading-edge  sweep  due  to  sideslip.  The  present  method  assumes  that,  for 
a given  angle  of  attack,  the  vortex  lift  remains  a constant  fraction  of 
the  attached  flow  lift  as  the  effective  leading-edge  sweep  changes. 

Hence,  for  fin  1,  which  has  an  increasing  effective  leading-edge  sweep 
as  increases,  the  theory  would  underpredict  the  data.  For  the  same 

reasons,  the  theory  would  overpredict  the  data  for  fin  3. 

The  discrepancy  between  data  and  theory  shown  in  figure  14  for 
fins  2 and  4 can  be  traced  to  measurement  error.  Fins  2 and  4 should 
have  the  same  lift  at  0=0.  If  the  data  of  fin  4 are  adjusted  so  that 
they  match  the  data  of  fin  2 at  0=0,  the  agreement  for  fin  4 is  good. 
The  agreement  for  fin  2 is  within  the  apparent  measurement  error. 

3.2.2  Aspect  ratio  1.0  fins.-  Data  for  the  Tj^^  set  of  fins  of 
body-tail  configuration  no.  1 were  obtained  for  several  Mach  numbers  and 
bank  angles  for  model  angles  of  attack  up  to  45°  ^ref.  16).  Since  the 
body  length  in  front  of  the  fins  was  large  enough  to  produce  strong  body 
vortices,  we  will  give  results  for  the  windward  fins  only.  The  fins 
were  not  deflected.  The  wing-alone  characteristics  for  angles  of  attack 
up  to  30°  were  obtained  by  interpolation  of  data  at  several  Mach  numbers 
for  fins  mounted  on  a reflection  plane  (ref.  17).  The  data  were  extra- 
polated to  higher  angles  of  attack  using  the  methods  of  reference  13. 

The  wing-alone  results  are  given  in  figure  15. 

Comparisons  of  theory  and  data  for  = 2.0  are  presented  in 
figure  16.  The  phenomena  discussed  above  for  the  1.3  aspect  ratio  fins 
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Figure  15.  Wing-alone  characteristics  for  low-aspect-ratio  fins  (^,4). 
” 2.0;  « 1.0,  Sjj  ” planform  area. 
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are  apparent  here  also.  Since  the  fins  have  delta  planforms,  the 

discussion  above  concerning  leading-edge  vortex  lift  is  applicable.  It 
is  interesting  to  note  that  for  * 30°  and  0 = 30°  fin  4 is  side- 
slipped 15°.  Since  the  fin  leading-edge  sweep  is  76°,  the  effective 
sweep  is  51°.  For  » 2.0,  it  seems  likely  that  the  flow  over  the 
leading  edge  is  supersonic.  Under  such  conditions  no  vortex  lift  could 
be  generated.  Thus,  the  correct  wing-alone  normal-force  curve  for  this 
case  would  be  for  the  attached  flow  only.  For  fin  4 at 

0 » 30°,  Ogq  = 45°.  If  we  extrapolate  the  slope  at  a^,  = 0 of  the  curve 
of  figure  15(a),  we  obtain  an  attached  flow  of  0.9  rather  than  the 

value  of  1.15  predicted  by  the  present  method.  The  improvement  in  the 
agreement  between  data  and  theory  indicates  that  the  above  explanation 
accounts  for  the  discrepancies  shown  in  figure  16.  It  thus  appears 
necessary  for  very  low-aspect-ratio  fins  to  consider  the  change  in 
leading-edge  and  side-edge  sweep  angles  in  calculating  vortex  lift. 


4.  COMPARISONS  WITH  DATA  FOR 
BODY-TAIL  CONFIGURATIONS 


4.1  Results  for  Bodies  Alone 

Comparisons  of  the  theory,  which  uses  equation  (1) , with  data  for 
bodies  alone  have  been  made  for  various  fineness  ratios  at  several  Mach 
numbers.  The  results  shown  in  figure  17  are  typical.  The  data  shown  are 
from  references  16  and  18.  Agreement  is  generally  good  except  when  the 
crossflow  Mach  number  is  transonic.  Either  the  full  two-dimensional 
crossflow  drag  is  not  applicable  for  yawed  bodies  at  transonic  speeds 
or  the  transonic  values  are  suspect.  Recent  experimental  results 

obtained  at  NASA/Ames  Research  Center  (ref.  19)  indicate  that  the 
values  are  reasonable.  Hence,  the  former  explanation  appears  most 
likely. 

In  the  present  method  the  body  vortex  strengths  and  positions  are 
based  on  a model  of  two  concentrated  vortices.  Recently,  however, 

Fidler,  Schwind,  and  Nielsen  (ref.  20)  have  shown  that  such  a model  is 
an  oversimplification  since  the  body  vorticity  can  be  distributed  over  a 
large  area.  The  reader  should  keep  this  in  mind  when  interpreting  the 
results  given  below  for  vortex-fin  interaction. 

4.2  Fin  Loads  in  the  Presence 
of  Body  Vortices 

The  purpose  of  this  section  is  to  uncover  phenomena  not  included  or 
not  properly  accounted  for  in  the  present  method  by  systematic  comparison 
of  data  and  theory.  We  will  be  especially  interested  in  phenomena  which 
become  more  important  at  higher  angles  of  attack  since  they  will  be 
pertinent  to  the  second  year's  work.  Results  for  each  of  the  leeward 
fins  and  for  the  windward  fins  are  presented  separately  below  to  facili- 
tate discussion  of  the  various  phenomena  of  interest.  Note  that  fin  1 rolls 
through  the  starboard  body  vortex  as  0 increases.  Fin  2 rolls  up  from  a 
horizontal  position.  It  is  never  "in" a vortex.  Fins  3 and  4 should  be 
weakly  affected  by  the  body  vorticity  for  positive  roll  angles  up  to  45°. 

4.2.1  Leeward  fins,  fin  1.-  In  addition  to  problems  involving  the 
use  of  a two-vortex  model  to  represent  distributed  body  vorticity,  the 
prediction  of  vortex-fin  interference  effects  must  take  into  account  the 
proximity  of  the  vortex  to  the  fin.  Figure  18  illustrates  this  phenomenon 
for  a rectangular  wing  near  a free  vortex  which  was  generated  upstream  by 


(d)  Vortex  far  from  wing;  small 
wing  thickness 


(a)  Vortex  far  from  wing;  large 
wing  thickness 


(e)  Vortex  close  to  wing;  small 
wing  thickness 


(b)  Vortex  close  to  wing;  large 


thickness 


(f)  Vortex  intersecting  wing 
small  wing  thickness 


(c)  Vortex  intersecting  wing 
large  wing  thickness 


18.  Illustration  of  vortex-wing  interference 
(photographs  taken  from  ref.  21). 


another  wing  (ref.  21)  . The  possibility  of  vortex  bursting  and/or  local 
wing  stall  is  handled  in  the  present  method  by  limiting  the  maximum  veloc- 
ity nearby  vortices  can  induce  on  the  fin.  A study  was  conducted  to 
determine  the  best  value  of  body-tail  configuration  no.  2 

shown  in  figure  3(d).  The  wing-alone  characteristics  for  the  fins 

(fig.  12)  were  used  since  the  T^  fins  closely  approximate  them  in 
planform.  The  results  of  the  study  for  fin  1, which  rolls  through  the 
right  body  vortex  as  0 increases  , are  shown  in  figure  19  for  = 1.75 
and  in  figure  20  for  •*  0.8  for  several  angles  of  attack.  The  two 
cases  show  the  effects  of  the  simple  expedient  of  limiting  the  velocity 
induced  on  the  fins  by  the  vortices.  The  best  values  of  V , /V„  appear 
to  be  0.4  for  M =1.75  and  0.1  for  M =0.8.  Since  V„_  /V  may  be 
thought  of  as  a local  stall  angle  (in  radians),  these  values  imply 
subsonic  stall  at  low  angles  of  attack  and  supersonic  stall  at  much 
higher  angles  of  attack. 

It  is  interesting  to  note  that  the  fin  1 normal  force  without  body 
vortices  present  is  underpredicted  for  both  = 1.75  and  = 0.8  for 
the  nearly  identical  fins  (see  figs.  13(a)  and  14(a)).  Hence,  we 

cannot  use  figure  19  and  20  to  determine  the  best  choise  of 
fin  1 results  for  rolling  moment  using  V _ /V„  = 0.4  for  M = 1.75  and 
0..1  for  = 0.8  are  presented  in  figures  21  and  22.  The  comparisons 
for  the  Mjjg  = 0,8  case  are  satisfactory.  However,  the  rolling  moments 
for  “ 1.75  are  not  well  predicted.  The  rolling  moment  together  with 
the  normal  force  of  a fin  in  a strong  encounter  with  a vortex  is  a sensi- 
tive indicator  of  the  performance  of  an  analytical  model  since  it  requires 
accurate  computations  of  the  center  of  pressure  location  as  well  as 
normal  force.  As  explained  in  reference  3 and  in  Appendix  B,  the  present 
theory  computes  the  location  of  the  center  of  pressure  using  reverse  flow 
theory  and  would  be  exact  if  the  flow  characteristics  were  linear.* 

Unknown  nonlinear  effects  not  accounted  for  by  reverse  flow  theory  may  be 
present.  Some  possible  effects  not  now  included  in  the  method  are 
(1)  vortex  bursting,  (2)  lateral  movement  of  a vortex  as  it  travels  over 
a fin,  (3)  increased  suction  due  to  the  greater  axial  velocity  in  the 
vortex  core,  and  (4)  nonlinear  terms  in  the  Bernoulli  pressure  equation. 


Note,  however,  that  the  present  code  assumes  that  the  interaction 
coefficients,  are  zero  when  i / j.  The  coefficient,  k.  is 

assumed  to  be  one.  See  Appendix  I', 
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19.  Effect  of  Vjjj^/Voo  on  prediction  of  \ 

fin  1 normal  force,  BT  2,  = 1.75;  J 

S„  = model  base  area.  J 
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Figure  20.  Effect  of  Vmav/V^  on  prediction  of 
fin  1 normal  force,  BT  2,  = 0.8; 

Sn  “ model  base  area. 
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Figure  22.  Comparison  of  data  and  theory  for  fin  1 
rolling  moment,  BT  2,  = 0.8,  = 0.1; 

Sj^  = model  base  area,  £ = model  base  diameter. 
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It  was  noted  that  the  measured  values  of  rolling  moment  shown  in 

figure  21  varied  with  0 in  a manner  similar  to  those  for  normal  force 

(fig.  19) . An  attempt  was  made  to  improve  the  rolling-moment  predictions 

for  = 1.75  by  using  the  wing-alone  center  of  pressure  location  given 

by  a_„.  The  results  are  shown  in  figure  23.  Note  that  the  agreement 
eq 

between  data  and  theory  improved  for  the  higher  angles  of  attack.  These 
results  seem  to  indicate  that  the  present  method  of  Appendix  B is  appro- 
priate if  the  vortex-fin  interaction  is  weak.  But,  if  the  interaction  is 
strong,  additional  factors  are  present  which  should  be  included  in  the 
theory . 

4.2.2  Leeward  fins,  fin  2.-  Comparisons  between  data  and  theory  for 

fin  2 of  body-tail  configuration  no.  2 which  rolls  up  from  a horizontal 

position  closer  to  the  port  body  vortex  are  given  in  figure  24  for 

M = 1.75  and  in  figure  25  for  M = 0.8.  See  sketch  2 below.  For  the 
00  -'00 


sketch  2.  Rear  view  of  body-tail  configuration  no.  2 
showing  position  of  body  vortices  in  crossflow 
plane  at  the  leading  edge  of  the  tail. 

supersonic  case,  the  predictions  are  somewhat  high  for  0 = 0°.  However, 
the  trends  are  well  predicted  including  the  fact  that  at  0 = 45°,  the 
normal  force  on  fin  2 decreases  with  increasing  body  angle  of  attack. 

For  the  svibsonic  case,  the  agreement  is  good  for  = 12.5°,  but 
for  the  higher  angles  of  attack  some  new  phenomenon  not  accounted  for  in 
the  theory  seems  to  be  present.  In  an  attempt  to  isolate  the  cause  of 
the  discrepancy,  the  fin  normal  force  was  calculated  for  the  case  of  no 
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body  vortices  present.  The  results  are  shown  in  figure  25  as  flagged 
symbols.  The  agreement  is  improved  somewhat  especially  at  0 = 45°. 
However,  for  the  smaller  values  of  <t>  and  for  = 16.6°  and  20.6°,  the 
normal  force  is  still  underpredicted  by  a significant  euxiount.  It  appears 
that  a source  of  favorable  lift  interference  is  present.  That  this  is 
so  can  be  seen  by  referring  to  figure  12.  (Note  that  the  reference  area 
for  figure  12  is  the  planform  area,  18.38  in.^,  and  for  figure  25  it  is 
the  model  base  area,  19.64  in.^.)  At  0 = 10°  in  figure  24,  the  meas- 
ured normal  force  for  = 20.6°  is  25  percent  higher  than  the 
wing-alone  curve  would  give  for  the  equivalent  angle  of  attack  indicated. 
It  appears,  then,  that  the  discrepancy  is  due  to  favorable  body  vortex- 
fin  interference  which  was  unsuspected  and  was  not  included  in  the  method. 

To  confirm  the  existence  of  the  phenomenon,  another  configuration 
with  a different  ratio  of  fin  semispan  to  body  radius  was  examined  also 
at  Mjjji  * 0.8.  The  configuration  is  body-tail  configuration  no.  3 shown 
in  figure  3(e).  The  wing-alone  normal  force  was  taken  from  reference  22 
zmd  is  presented  in  figure  26.  The  data  were  taken  from  reference  16  and 
are  presented  for  several  auigles  of  attack  in  figure  27.  Note  the 
similarity  between  these  data  and  the  data  of  figure  25.  The  values  of 
fin  normal  force  which  would  be  predicted  if  no  body  vortices  were  pre- 
sent are  also  shown  in  figure  27.  It  is  clear  that  favorable  vortex 
interference  is  present  for  0 ^ 0 < 20°.  It  is  interesting  to  note  that 
for  both  cases  the  magnitude  of  the  favorable  interaction  seems  to  be  a 
function  of  the  distance  of  the  body  vortex  from  the  fin.  In  fact,  when 
the  f'n  is  at  20°  bank  (second  quadrant) , the  favorable  interaction  is 
lost.  The  mechanism  here  appears  to  be  similar  to  that  encountered  for 
strake-wing  interference.  The  phenomenon  will  be  investigated  in  the 
second  year  of  the  present  effort. 

4.2.3  Windward  fins.-  Comparisons  between  data  and  theory  for 
fins  3 and  4 of  body-tail  configuration  no.  2 at  «*  0.8  are  given  in 
figure  28.  Note  that  the  trends  are  similrr  to  those  of  figure  14  as 
they  should  be  since  the  fins  are  nearly  identical  in  planform.  The 
chief  difference  is  the  size  of  the  fins  relative  to  the  body.  In 
section  3.2,  the  possibility  was  discussed  of  improving  the  prediction 
for  fin  3 by  taking  into  account  the  dependence  of  vortex  lift  on 
leading-edge  sweep  changes  due  to  pitch  and  roll. 
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Fin  angle  of  attack,  a,  degrees 


Figure  26.  Wing-alone  normal-force  coefficient  for 
low-aspect-ratio  fins  b 0.8; 

» - 1.0,  Sd  ” planform  area. 
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Figure  28.  Comparison  of  data  and  theory  for  fin  3 
normal  force,  BT  2,  = 0.8,  = 0.1; 

Sj^  = model  base  area. 


T^ese  improvements  are  based  on  an  approximate  method  now  described. 
For  clarity  the  discussion  will  be  made  initially  for  triangular  wings. 
The  normal-force  curve  for  a wing  alone  which  develops  vortex  lift  can 
be  approximated  by  the  following  formula. 


sin  a cos  a + 


0=0 


K Sin 

V 


o 


(34) 


The  corresponding  normal-force  curve  has  the  following  characteristic 
shape  : 


Sketch  3.  Normal-force  curve  for  a wing  alone 
which  develops  vortex  lift. 


where 


do  / „ 

/ opO 


normal-force  curve  slope  at  a 


0 


= "vortex  lift"  constant 


The  constant  can  be  obtained  from  an  experimental  normal-force  curve. 

It  can  also  be  estimated  theoretically  by  assuming  that  all  the  leading- 
edge  suction  IS  converted  into  vortex  lift  (Polhamus  analogy) . For  delta 
wings  at  subsonic  speeds  theoretical  values  of  depend  only  on 

leading-edge  sweep  angle  or  its  equivalent,  aspect  ratio.  For  supersonic 
speeds,  the  theoretical  value  of  depends  on  the  ratio  of  the  tangent 

of  the  leading-edge  sweep  angle  to  the  tangent  of  the  Mach  angle,  or 
3®.  For  supersonic  edges,  is  zero.  Generally,  not  all  the  leading- 
edge  suction  is  converted  to  vortex  lift.  The  parameter  K*  is  defined 
to  represent  this  fraction. 
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(35) 


K*  = exp 

theory 

Correlations  of  K*  for  delta  wings  are  included  in  reference  12. 

For  a cruciform  missile  at  included  angle  of  attack,  a^,  and  roll 
angle,  0,  the  sweep  angle  of  a delta  wing  is  changed  from  its  nominal 
angle  at  = 0.  This  changes  both  the  potential  normal  force  and  the 
vortex  normal  force.  The  change  in  normal  force  is  handled  in  the 
present  method  by  the  change  in  equivalent  angle  of  attack  determined 
using  the  interference  factor,  K^.  This  is  correct  for  the  normal  force 
due  to  potential  flow.  However,  special  consideration  must  be  given  to 
the  vortex  normal  force  because  not  only  is  the  leading-edge  suction 
changed  by  change  in  sideslip  angle  but  so  also  is  the  fraction  of  it 
converted  to  vortex  lift. 

Consider  first  the  changes  in  sweep  angle  of  a delta  fin  leading 
edge  as  influenced  by  pitch  and  roll.  For  zero  pitch  and  roll,  we  have 

Ao  = cot-"  (36) 

For  combined  pitch  and  roll,  the  sweep  angle,  A^^^,  is  now 

^eff  = (®)-  ^ 

From  the  results  of  reference  12,  values  of  K*  corresponding  to  A^ 
and  can  be  obtained  and  are  designated  K*  and  K*^^. 

In  estimating  the  change  in  vortex  normal  force,  we  determine 
in  the  usual  way.  The  curve,  (dCjj^/da)  ^Qsin  cos  a,  is  constructed  on 
the  normal-force  curve  of  the  wing  alone  as  shown  in  the  preceding  sketch. 
A "vortex  normal  force"  is  determined  by  subtracting  the  potential  nor- 
mal force  from  the  experimental  normal  force.  This  vortex  normal  force 
is  then  scaled  by  the  ratio,  This  method  assumes  that  the 

leading-edge  suction  is  a constant  fraction  of  the  potential  normal 
force  as  sweep  angle  changes.  Using  vortex-lattice  or  panel  methods, 
we  can  determine  the  actual  change  in  leading-edge  suction  with  sweep 
angle  before  applying  the  K*  factors  but  have  not  done  so  in  the 
approximate  method. 
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An  application  of  the  above  improvement  in  the  present  method  has 
been  made  to  fin  1 for  the  conditions  of  figure  28(a) . The  fin  is  a 
clipped  delta  fin,  and  it  has  been  approximated  by  a delta  fin  with  the 
same  aspect  ratio  for  obtaining  values  of  K*.  In  determining  the  poten- 
tial normal  force,  sin  a was  taken  as  a and  cos  a as  unity  in  equa- 
tion (34).  The  pertinent  wing  normal-force  curve  is  given  in  figure  12(a). 
The  numerical  results  obtained  by  this  "improved  method"  are  shown  in 
figure  28(a)  by  the  flagged  symbols.  Much  improved  agreement  with  data 
is  exhibited.  It  seems  clear  from  the  improvement  shown  in  figure  28  (a) 
that  an  important  phenomenon  not  now  included  in  the  present  method  has 
been  identified.  Furthermore,  a means  for  accounting  for  its  effects 
has  been  found.  Further  refinement  to  account  separately  for  the 
effects  of  leading  edges  and  side  edges  seems  desirable. 

4.3  Total  Loads  for  Body-Tail  Configurations 

We  have  made  comparisons  of  data  for  complete  body-tail  configurations 
with  the  present  theory  without  the  improvements  discussed  above.  The 
results  shown  in  figures  29-31  are  typical.  We  will  discuss  the  normal- 
force  results  first  and  then  rolling  moment  and  side  force. 

Generally,  the  agreement  for  total  normal  force  is  good  at  low 
and  fair  at  high  a^.  Let  us  first  consider  the  results  for  body-tail 
configuration  no.  2 at  = 1.75.  Note  that  the  contribution  of  the 
tail  fin  loads  to  the  normal  force  is  well  predicted.  Results  not  pre- 
sented show  that  the  body-alone  normal  force  for  this  case  is  also  well 
predicted.  Consideration  of  the  axial  location  of  the  center  of  pressure 
(fig.  29(b))  also  suggests  that  the  predicted  contribution  to  the  normal 
force  for  the  tail  region  at  0 = 45°  is  slightly  low  since  the 
predicted  center  pressure  location  is  slightly  forward. 

For  the  M = O.o  case,  it  is  clear  from  fiqvr.a  30(a)  that  the 
calculated  contribution  of  the  tail  fins  to  normal  force  is  somewhat  low 
at  high  angles  of  attack.  As  discussed  above  in  section  4.2.2,  the 
discrepancy  appears  to  be  due  to  favorable  vortex-fin  interaction  which 
is  not  accounted  for  in  the  present  method. 

Results  for  the  AIM-9L  body  and  tail  (body-tail  configuration  no.  4) 
are  shown  in  figure  31  for  = 0.8  and  a bank  angle  of  22.5°.  The 
fins  were  undeflected.  The  agreement  between  data  and  theory  for  normal 
force  is  very  good. 
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The  induced  rolling  moments  for  the  body-tail  cases  presented  in 
figures  29(c)  emd  30(c)  and  31(b)  are  interesting  since  the  rolling 
moment  is  a second-order  quantity  (linear  theory  with  no  vortices  would 
give  zero  rolling  moment) . The  rolling  moments  shown  correspond  roughly 
to  2°  of  aileron  control.  The  agreement  between  data  and  theory  is  only 
fair.  The  agreement  should  improve  considerably  if  the  phenomena 
discussed  in  sections  3 and  4.2  are  incorporated  into  the  method. 
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Included  angle  of  attack,  a^,  degrees 
(a)  Normal-force  coefficient  in  unrolled  body  coordinates 


Figure  29. 


Comparison  of  data  and  theory  for  BT  2,  M 
V^ax/'^oc  “ 0.4;  Sr  = model  base  area,  “ 
= model  base  diameter. 


Included  angle  of  attack,  a^,  degrees 
(c)  Rolling-moment  coefficient. 
Figure  29.  Concluded. 
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(a)  Normal-force  coefficient  in  unrolled  body  coordinates. 


Figure  30.  Comparison  of  data  and  theory  for  BT  2,  = 0.8, 

= 0.1;  Sr  = model  base  area, 

= model  base  diameter. 
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Included  angle  of  attack,  a^,  degrees 

(b)  Rolling-moment  coefficient. 

Figure  31.  Comparison  of  data  and  theory  for  BT  4, 
= 0.8,  Vniax/^oc.  = 0.4;  = model  base  area, 

£ = model  base  diameter. 


5.  CCMPARISONS  WITH  DATA  FOR  BODY -CANARD  AND 
BODY-CANARD-TAIL  CONFIGURATIONS 

This  section  is  concerned  with  body-canard  and  body-canard-tail 
configurations  which  require  vortex  tracking  for  the  calculation  of 
forces  and  moments.  The  method  used  to  obtain  the  vortex  system  behind 
the  canard  fins  is  heuristic,  and  its  validity  can  best  be  determined  by 
comparing  its  predictions  with  data.  The  method,  as  described  in 
Appendix  D,  represents  any  body -generated  vorticity  which  may  form  aft 
of  the  canard  trailing  edge  as  discrete  vortices  which  change  strength 
and  position  as  a function  of  their  axial  positions.  The  first  sub- 
section below  contains  results  for  body-canard  configurations  (see  fig.  3). 
Since  no  tail  fins  are  present,  these  results  are  a sensitive  indication 
of  the  ability  of  the  tracking  method  to  reproduce  actual  body  loads  in 
the  presence  of  the  weike  over  the  afterbody.  The  second  subsection  below 
is  concerned  with  body-canard-tail  configurations  (see  fig.  3) . These 
results  are  good  tests  of  the  ability  of  the  track..ng  method  to  locate 
the  vortices  of  the  afterbody  wake  properly  with  respect  to  the  tail  fins. 

5.1  Results  for  Body -Canard  Configurations 

One  of  the  features  of  the  present  method  is  that  all  vorticity  is 
regarded  as  discrete  and  all  the  vorticity  generated  from  a separation 
line  on  the  body  is  represented  by  a single  vortex.  For  example,  the 
nose  vorticity  0’~  the  port  side,  the  trailing-edge  vorticity  from  canard 
fin  1 and  the  afterbody  vorticity  on  the  port  side  are  each  represented 
by  a single  vortex.  Vapor  screen  (ref.  5)  and  laser  anemometer  (ref.  20) 
studies  of  the  development  of  such  vorticity  indicate  that  more  than  one 
vortex  may  at  times  be  req'uired  to  represent  the  shed  vorticity  from  each 
of  these  sources.  It  is  known,  for  example,  from  the  work  of  reference  5 
that  for  < 20°  the  vorticity  formed  ahead  of  the  canard  fins  of 
body -canard  configuration  no.  1 is  dispersed  as  it  moves  over  the  fins 
and  is  no  longer  rolled  up  by  the  time  it  reaches  the  canard  trailing 
edge.  The  present  model  represents  such  a vorticity  field  by  a concen- 
trated vortex  which  cannot  disperse.  The  user  of  the  code  does  have  the 
option,  however,  of  eliminating  the  nose  vortices  at  the  canard  leading 
or  trailing  edges  if  he  feels  that  they  have  dispersed  in  the  actual  flow. 
The  vapor  screen  study  of  reference  5 also  implied  that  vorticity  shed 
from  one  of  the  windward  fins  can  travel  up  and  around  the  body  and 
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strongly  interact  with  the  body  boundary  layer.  The  method  used  herein 
does  not  include  this  effect. 

During  the  tests  of  reference  4,  individual  fin  load  data  were 
obtained  for  body-canard  configuration  nos.  1 and  2 shown  in  figure  3. 
These  data  enable  us  to  investigate  the  effects  of  afterbody  wake  asym- 
metry on  total  body  loads.  Results  for  body-canard  configuration  no.  1 
are  presented  in  figure  32  for  = 1.75  and  zero  bank.  The  moderate- 
aspect-ratio  Cg  fins  are  deflected  15°  in  yaw  and  0°  in  pitch.  The 
comparisons  show  that  total  normal  force  is  accurately  predicted. 

However,  side  force  is  poorly  predicted  although  the  trend  is  correct. 
Note  that  the  contributions  of  the  fins  to  the  total  normal  and  side 
forces  are  well  predicted.  Hence,  any  discrepancies  present  must  be 
attributed  to  the  afterbody  wake  model. 

Results  for  body-canard  configuration  no.  2 are  presented  in  fig- 
ure 33  for  Mjjjj  = 1.75  and  20°  bank.  The  low-aspect-ratio  fins 

were  undeflected.  Note  that  the  Cy  fins  are  considerably  smaller 
relative  to  the  body  than  the  Cg  fins,  yet  similar  results  were 
obtained.  Again,  the  contribution  of  the  fins  to  total  loads  is  well 
predicted.  The  agreement  for  total  normal  force  is  satisfactory  though 
not  quite  so  good  as  for  the  previous  case.  The  results  for  total  side 
force  are  interesting  since,  as  in  the  previous  case,  the  trend  (for 
< 24°)  is  represented  correctly  although  the  magnitude  is  not.  In 
this  case,  the  side  force  is  a second-order  quantity  being  generally  less 
than  10  percent  of  the  normal  force. 

Another  set  of  comparisons  are  presented  in  figure  34  for  body- 
canard  configuration  no.  3 for  = 0.8  and  zero  bank.  This  configu- 
ration represents  the  AIM-9L.  Note  that  the  canards  (see  fig.  3(g))  are 
considerably  larger  relative  to  the  body  than  the  fins  of  body-canard 
configuration  no.  1.  Note  also  that  the  leading  edge  of  each  canard  fin 
has  a break  in  sweep.  Although  the  present  version  of  the  code  cannot 
handle  breaks  in  leading-edge  sweep,  the  fins  can  be  adequately  modeled 
with  a straight  leading  edge  running  from  the  actual  location  of  the 
leading  edge  of  the  root  chord  to  the  actual  location  of  the  tip  chord. 
For  the  results  of  figure  34,  the  fins  were  deflected  20°  in  yaw  and  0° 
in  pitch.  The  agreement  is  fair  for  normal  force  and  poor  for  side 
force.  (Note  that  the  side-force  scale  is  50  times  the  normal-force 
scale.)  If  more  accurate  side-force  calculations  on  a percentage  basis 
are  to  be  made,  some  improvement  in  the  afterbody  vortex  model  is 
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Included  angle  of  attack,  a^,  degrees 
(a)  Normal-force  coefficient  in  unrolled  body  coordinates 


Figure  32.  Comparison  of  data  and  theory  for  BC  1,  - 1.75 

V /V  - 0.4;  S„  - model  base  area. 
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Included  angle  of  attack,  a^,  degrees 
(b)  Side-force  coefficient  in  unrolled  body  coordinates 


Figure  32.  Concluded. 


indicated.  The  sources  of  side  force  can  be  found  by  examining  the 
vortex  impulse  theorem  (eq.  (28))  which  is  used  to  compute  them.  If 
we  tckke  the  real  and  imaginary  parts  of  equation  (28)  and  integrate  over 
the  afterbody  section,  we  obtain 

, NV(x)r r / V n 


body 

base 


L \ v,3  /J 


trailing  edge  of 
canard  section  ' 


NV(x) 


' 4.3) 


body 

base 


trailing  edge  of 
canard  section 


where  F ' is  nondimensionalized  by  2vVa  and  y and  z are  nondimen- 
sionalized  by  a.  If  a vortex  is  free,  its  strength  is  constant,  and  its 
contribution  to  the  normal  force  and  side  force  depend  on  its  movement 
in  the  y^  and  directions,  respectively.  The  initial  strength  of  an 
afterbody  vortex  is  zero  and  its  contribution  to  the  forces  depends  only 
on  its  strength  and  position  in  the  crossflow  plane  at  the  body  base. 

We  have  computed  the  predicted  effect  of  each  vortex  present  on  the 
afterbody  loads  for  the  case  shown  in  figure  32  at  = 20.6°.  The 
results  are  given  in  Table  1.  Note  that  a 5-  to  10-percent  change  in 
the  vertical  position  relative  to  the  body  axis  of  any  of  the  vortices 
could  make  a large  difference  in  the  predicted  side  force  for  the  after- 
body section.  Also  large  side-force  contributions  of  alternating  sign 
occur.  The  results  of  Table  1 clarify  the  need  for  accurate  vortex 
tracking  past  the  afterbody  section.  While  the  present  mode'^  is  capable 
of  reproducing  most  of  the  physics,  further  improvement  is  needed  if 
accurate  results  for  side  force  are  to  be  obtained. 


5.2  Results  for  Body-Canard-Tail  Configurations 

Comparisons  between  data  emd  theory  for  two  different  body-canard- 

tail  configurations  are  given  in  this  section.  The  configurations  are 

unrolled  but  have  large  yaw  control  deflections.  Figure  35  presents 

results  for  body-canard-tail  configuration  no.  1 at  M - 1.75  with 

00 

zero  bank.  The  fins  are  deflected  15°  in  yaw  and  0°  in  pitch.  The 
normal  force  is  well  predicted.  However,  the  theory  for  side  force 
departs  increasingly  from  data  for  angles  of  attack  above  18°.  Since 
the  contribution  of  the  fins  to  the  total  side  force  is  fairly  well 
predicted,  the  procedure  for  obtaining  the  afterbody  vortex  positions 
and  strengths  on  the  afterbody  section  is  the  likely  sovirce  of  the 
discrepancy.  The  results  for  rolling  moment  are  only  fair. 

Since  the  body  cannot  support  a rolling  moment,  an  analysis  of  the 
individual  contributions  of  the  fins  to  the  total  rolling  moment  should 
be  of  interest.  The  results  are  given  in  Table  2 for  a^,  “ 8.5°  and  20.6°. 
It  is  interesting  to  note  that  both  the  canard  section  and  the  tail 
section  contribute  substemtially  to  the  rolling  moment  but  for  different 
reasons.  The  canard  section  does  not  have  significeint  body  vorticity 
present.  Hence,  the  only  source  of  rolling  moment  for  the  canard  section 
is  differential  fin  loading  due  to  sideslip.  Canard  fins  2 and  4 are  not 
sideslipping  as  increases.  Therefore,  fins  1 and  3 contribute  most 

to  the  total  rolling  moment. 

In  the  absence  of  vorticity  shed  by  the  body  and  canard  fins,  the 
tail  section  for  0=0  cannot  support  a rolling  moment.  Hence,  the 
only  source  of  rolling  moment  for  the  tail  section  is  vortex-fin  inter- 
ference. This  is  clearly  evident  in  the  data  for  fins  1,  2,  and  4.  Note 
that  the  discrepancies  between  data  and  theory  for  individual  fins  are  not 
large  compared  to  the  size  of  the  loads  seen.  Yet  those  same  discrepancies 
create  large  percentage  errors  in  the  overall  rolling  moment. 

It  is  interesting  to  see  which  of  the  vortices  present  in  the  wake 
influence  each  fin.  Calculated  and  measured  vortex  positions  for  several 
angles  of  attack  at  the  leading  edge  of  the  tail  root  chord  are  presented 
in  figure  36.  Note  that  the  predictdons  of  the  present  method  place  the. 
left  afterbody  vortex  too  far  to  the  right  and  the  right  afterbody  vortex 
too  high. 
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TABLE  2.  ANALYSIS  OF  ERROR  IN  ROLLING  MOMENT  FOR 

BODY-CANARD-TAIL  CONFIGURATION  NO.  1‘, 

M^  = 1.75,  0 = 0,  6^  = 63  = 15°, 

6„  = 6^  = 0°,  S = model  base  area, 

24  R 

- model  base  diameter. 

(a)  = 8.5°. 


~ data 

theory  ~ ^^i^data]^^^  ^ 


Data  Theory  < 


all  fins 


1 CRMCl 

0.1605 

0.1634 

0.0029 

0.06 

CRMC2 

0.1423 

0.1504 

0.0081 

0.17 

CRMC3 

0.1974 

0.1829 

-0.0145 

-0.30 

CRMC4 

0.1392 

0.1545 

0.C153 

0.32 

CRMTl 

-0.0747 

-0.0575 

CM 

0 

• 

0 

0.36 

CRMT2 

0.3170 

0.3454 

0.0284 

0.59 

CRMT3 

-0.0001 

0.0118 

0.0119 

0.25 

CRMT4 

0.2457 

0.2924 

0..0467 

0.97 

canard 

0.0400 

0.0154 

fins 

^^i^tail 

0.1459 

0.1223 

fins 

0.1859 

0.1377 

fins 
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theory 

“"'’■'“•I  canard 
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theory 
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(b)  = 20.6°. 


Data 

Theory 

^i 

CRMCl 

0.1512 

0.1494 

-0.0018 

-0.01 

CRMC2 

0.2775 

0.3018 

0.0243 

0.16 

CRMC3 

0.2395 

0.1968 

■0.0427 

-0.29 

CRMC4 

0.2854 

0.3102 

0.0248 

0.17 

CRMTl 

0.1446 

0.1272 

-0.0174 

-0.12 

CRMT2 

0.6094 

0.7496 

0.1402 

0.95 

CRMT3 

-G.0044 

0.043^ 

0.0478 

0.32 

CRMT4 

0.6083 

0.6250 

0.0167 

0.11 

<=<>ca„ard 

0.0391 

fins 

-0.1479 

0.0408 

fins 

-0.0675 

0.0799 

fins 

theory  ~ 

^^‘^^t^]canard 

fins 

theory 
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fins 
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In  the  present  method,  the  afterbody  vortices  are  assumed  to  be 
initially  located  50°  on  either  side  of  the  crossflow  velocity  vector 
given  by  the  free-streeun  flow  and  by  the  vorticity  shed  by  the  canard 
section.  For  the  case  presented  in  figure  36,  the  crossflow  velocity 
vector  at  the  trailing  edge  of  the  canard  fins  points  to  the  right  25°, 
19°,  and  15°  from  the  vertical  for  a^.  ■ 12.5°,  16.6°,  and  20.6°, 
respectively.  Ihis  means  that  the  left  afterbody  vortex  would  be 
started  in  the  second  quadrant  at  65°,  59°,  and  55°  measured  clockwise 
from  fin  T2.  However,  the  vapor  screen  results  of  reference  5*  indicate 
that  the  left  afterbody  vorticity  actually  separates  from  the  body  much 
lower  in  the  second  quadrant.  Furthermore,  the  angular  location  of  the 
separation  line  does  not  vary  signif icemtly  with  a^. 

An  attempt  was  made  to  represent  more  accurately  the  physics  for 
this  case  by  assuming  the  afterbody  vorticity  to  separate  60°  on  either 
side  of  the  vertical.  The  results  for  the  modified  theory  are  presented 
in  figure  36  as  flagged  symbols.  Note  that  the  agreement  between  data 
and  theory  for  is  considerably  better.  No  improvement  was  obtained 

for  A„,  however.  The  data  of  reference  5 showed  that  the  C3  vortex 
moved  up  and  to  the  right  as  it  moved  downstream.  At  some  point  along 
the  body  depending  upon  a^,  it  came  close  enough  to  the  body  surface 
to  interact  with  the  body  boundary  layer.  From  that  point  aft,  it  was 
combined  with  the  right  afterbody  vortex.  That  behavior  is  not 
accounted  for  in  the  present  method. 

In  order  to  show  that  the  results  discussed  above  were  not  limited 
to  an  isolated  case,  comparisons  between  data  and  theory  were  made  for 
another  configuration  and  are  presented  in  figure  37.  The  model  is  the 
AIM*9t.  of  reference  14  (body-canard-tail  configuration  no.  2,  fig.  3(i)). 
The  model  is  unrolled  and  has  20°  of  yaw  control  and  0°  of  pitch  control. 
The  free-3tre2un  Mach  number  is  0.8.  As  in  the  case  previously  discussed, 
the  agreement  is  good  for  normal  force  but  only  fair  for  side  force  and 
rolling  moment. 


See  also  figure  20  of  reference  4. 


6.  CONCLUDING  REMARKS 


This  report  covers  the  first  year's  work  in  a two-year  effort  to 
develop  a computer  progreun  for  determining  the  forces  and  moments  on  the 
individual  components  of  cruciform  missiles  at  hi-^h  angles  of  attack. 

The  program  is  intended  for  use  at  subsonic,  transonic,  and  supersonic 
speeds.  It  presently  accommodates  arbitrary  roll  angles,  pitch  control, 
and  yaw  control.  It  is,  however,  subject  to  the  following  limitations. 

1.  The  fins  in  each  finned  section  must  be  identical,  unceuiibered, 
and  untwisted. 

2.  Opposing  fins  must  be  symmetrically  deflected;  the  program  does 
not  accommodate  roll  control  at  this  time. 

3.  The  wake  from  each  fin  is  represented  by  a single  vortex. 

4.  The  trailing  set  of  fins  may  not  be  deflected. 

5.  The  leading  edges  of  the  fins  must  not  be  sweptf orward , and  the 
trailing  edges  must  be  unswept. 

6.  The  small  angle  definitions  of  a and  p are  used. 

7.  The  body  wake  upstream  of  the  first  finned  section  is  represented 
by  two  discrete  symmetric  vortices. 

8.  The  vorticity  shed  from  the  afterbody  section  is  represented  by 
two  discrete  vortices. 

As  part  of  this  effort,  a comprehensive  comparison  between  data  and 
predictions  of  the  computer  program  has  been  made  for  bodies  alone,  body- 
tail  combinations,  body-canard  combinations,  and  body-canard-tail 
combinations.  The  following  ranges  of  parameters  were  covered: 


Angle  of  Attack: 

0°  - 24° 
0°  - 47° 

for 

for 

finned 

bodies 

missiles 

alone 

Mach  Number: 

0.8  - 2.0 

Roll  Angle: 

o 

in 

1 

o 

o 

Control  Deflecv.ion: 

0°  - 20° 

Fin  Aspect  Ra'.io: 

0.87  - 4.2 
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Hie  following  conclusions  can  be  drawn  from  the  study; 

1.  For  the  bodies  alone,  accurate  results  were  obtained  except 
when  the  crossflow  Mach  nuirber  was  near  unity.  The  two-dimensional 
crossflow  drag  coefficients  used  in  the  procedure  are  too  large  under 
these  circumstances. 

2.  For  computing  individual  fin  loads,  the  equivalent  cuigle  cf 
attack  concept  is  generally  adequate  except  for  two  new  vortical  effects 
not  included  in  the  method  as  described  in  the  next  two  items. 

3.  As  a missile  rolls,  its  fin  leading  and  trailing  edges  undergo 
sweep  changes.  Ihese  sweep  changes  cause  changes  in  the  leading-edge 
and  side-edge  suctions  as  well  as  the  fractions  of  these  suctions 
converted  to  vortex  lift.  This  effect,  not  included  in  the  computer 
program,  can  cause  significant  changes  in  panel  normal  force  particularly 
for  low-aspect-ratio  fins. 

4.  Increased  panel  normal  force  can  result  under  certain  conditions 
as  a result  of  possible  favorable  vortex  interference.  This  occurred  for 
the  left  horizontal  fin  (T2)  in  the  16°  to  29°  range  for  angles  of  roll 
between  0°  and  30°  with  low-aspect-ratio  fins. 

5.  In  general,  normal  force  and  axial  location  of  center  of  pressure 
were  well  predicted  by  the  method  for  complete  configurations. 

6.  For  no  control  deflections,  a cruciform  configuration  can  develop 
side  forces  and  rolling  moments  principally  as  a result  of  nose,  after- 
body, or  canard  vortices.  In  the  rsinge  of  cuigle  of  attack  considered 
here,  the  side  force  is  typically  about  10  percent  of  the  normal  force, 
and  the  rolling  moment  corresponds  to  2°  of  roll  control  on  all  four 
fins.  The  errors  in  predicting  these  quantities  are  of  the  order  of 

50  percent.  Improving  the  accuracy  of  these  predictions  will  depend  to 
a considerable  degree  on  developing  a better  model  for  predicting  the 
vortex  path  over  the  afterbody  and  the  afterbody  vortex  strengths  and 
predictions. 

In  the  second  year's  work  we  plan  to  compare  the  theory  with  addi- 
tional diagnostic  data,  particularly  that  which  will  be  obtained  to 
angles  of  attack  of  about  50°  in  the  Ames  Research  Center  11-Foot  Wind 
Tunnel  under  a MASA/Ames  Research  Center  contract.  Changes  in  the 
computer  program  to  improve  its  accuracy  at  high  angles  of  attack  are 
contemplated,  euid  certain  of  the  limitations  and  shortcomings  discussed 
above  will  be  removed. 


Ill 
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APPENDIX  A 

PROCEDURE  FOR  COMPUTATION  OF  VORTEX 
PATHS  IN  OHE  PRESENCE  OF  A CRUCIFORM 
WING-BODY  COMBINATION 


The  determination  of  the  loading  on  a finned  section  of  a missile 
due  to  the  presence  of  a set  of  vortices  requires  that  the  paths  of  the 
vortices  be  known.  A procedure  for  computing  those  paths  has  been  pre- 
viously derived  (refs.  9,  12,  and  A.l)  and  is  repeated  below.  It  is 
based  on  slender-body  theory  and  assumes  that  the  vortices  can  be 
"tracked"  in  a crossflow  plane  which  is  perpendicular  to  the  body  axis 
and  which  moves  downstre2un  at  V^cos  a^.  The  body  boundary  conditions 
are  properly  handled  by  imaging  the  external  vortices.  Since  the  image 
system  for  a cruciform  wing-body  combination  would  be  quite  complex,  the 
actual  configuration  in  the  crossflow  plane  is  transformed  into  a circle 
as  shown  in  the  sketch  below.  The  sketch  shows  two  asymmetrical  vortices 
in  the  presence  of  a cruciform  wing-body  combination  at  combined  pitch 
and  roll. 


(a)  Physical  plzme. 


(b)  Tramsform  plane. 


Sketch  A.l  Computational  planes  for 
vortex  tracking. 


l^e  vortices  are  assumed  to  be  discrete  potential  vortices  moving  as 
Lagrangian  fluid  particles.  We  will  derive  the  equation  of  motion  of 
Fj  with  Fg  present,  and  then  generalize  the  results  to  many  vortices. 
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The  fin  semispan,  s,  will  be  considered  to  be  a fiaiction  of  the  distance 
along  the  body  axis,  x.  No  effects  of  fin  deflection  on  vortex  paths 
are  included. 


Consider  first  the  transformation  and  its  derivatives  of  the  cruci 
form  wing -body  combination  in  the  a plane  into  the  circle  of  radius, 
r in  the  v plane. 


The  complex  potential  for  the  flow  in  the  v plane  is  given  by 


We  have  assvuned  that  the  body  is  of  uniform  radius  and  have  put  no  image 
vortices  at  the  center  of  the  body.  The  direction  of  the  flow  at 
infinity  is  unchanged  under  the  transformation. 

The  complex  potential  W^^fvCa)]  for  the  complex  velocity  v^^  - iw^ 
of  vortex,  , in  the  a plane  is  the  total  complex  potential  W(v) 
for  the  flow  including  the  vortices  minus  the  singularity  at  vortex,  F . 


r 1 fo  .4  V ( 

3 

1 - aVaV. 

If  there  are  n external  vortices,  , . . . , r^,  then  equation  (A. 10) 

will  have  the  form 


V - xw 

1 1 


(A. 11) 


For  the  j vortex,  we  can  write 


. -uL  . 1 \ 


V . - iw  . = -iV„sin 
ID 


■ V k / 1 

^ Z_.  2tt  I V . - V, 

v=i  \ J ^ 


1 /da\/d^v  \ 

2 /-  / \ day„  „ 27t  \dv  / V , 2 / 


(A. 12) 

In  determining  the  trajectories  of  vortices  passing  by  a cr'ciform 
finned  section  of  a missile,  the  positions  and  strengths  of  the  vortices 
must  be  given  at  some  initial  value  of  x.  Then  the  appropriate  equations 
of  motion  for  the  vortices  are 


dy^  dz . 

= V . and  = w . 

dt  j dt  3 


(A. 13) 


Since  the  crossflow  plane  is  moving  downstream  with  the  speed,  V^cos  a^, 
we  have 


dx 

^ = V cos  a 
dt  00  c 


(A.  14) 


Substituting  equation  (A  .14)  into  (A  .13)  gives 


^ Xi  1 
_ ■■■  « ■ - 
dx  V„  cos 


^ = !!i  1 

dx  V„  cos  a 
00 


(A. 15) 


120 


The  set  of  differential  equations  (A. 15)  together  with  the  initial  condi- 
tions and  the  set  of  algebraic  relations  (A. 12)  form  a system  of  ordinary 
first-order  initial  value  differential  equations  which  can  be  solved  by 
standard  integration  routines. 

The  present  method  uses  small  angle  approximations  throughout. 

Hence,  the  quantity  cos  in  equations  (A. 15)  is  set  equal  to  one. 

When  the  influence  of  the  fins  is  to  be  ignored,  the  code  sets  s(x)  = a. 

The  following  relationships  can  be  used  to  go  back  and  forth  between 
the  a and  v planes: 


(v^  + r^v^)  + V(v^  + -” 


(a^  + a*/a^)  + ~\}  + t-Vcr")^  - 4r^ 




(A.  16) 


(A. 17) 


The  roots  are  to  be  chosen  such  that  corresponding  points  are  always  in 
corresponding  quadrants. 
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APPENDIX  B 


REVERSE  FLOW  METHOD  FOR  COMPUTATION 
OF  FIN  FORCES  AND  MOMENTS  IN 
THE  PRESENCE  OF  VORTICES 

B . 1 . INTRODUCTION 

The  purpose  of  this  appendix  is  to  present  the  procedure  used  in 
the  cede  MLOADS  to  compute  the  effects  of  vortices  on  the  loading  exper- 
ienced by  a cruciform  wing-body  combination.  The  reverse  flow  method 
used  herein  was  originally  developed  for  the  calculation  of  forces  and 
moments  on  an  external  store  in  a nonuniforra  flow  field  and  was  presented 
in  reference  B.l.  The  development  below  closely  follows  that  of 
Appendix  I of  reference  B.l. 

In  the  procedure  under  discussion  the  body  is  at  zero  angle  of 
attack  and  the  fins  are  undeflected  so  that  only  the  loading  due  to 
vortices  will  be  computed. 


B.2.  DISCUSSION  OF  METHOD 


Our  problem  is  to  determine  the  fin  normal  forces  and  rolling 
mome^'ts  due  to  an  external  vortex  passing  a cruciform  fin-body  combina- 
tion. The  first  step  is  to  determine  the  changes  in  fin  twists  due  to 
the  vortex.  The  vortex  is  assumed  to  be  rectilinear  and  parallel  to 
the  fin,  and  it  induces  an  angle  of  attack  normal  to  the  i'th  fin, 
a . (t)  which  can  be  established  by  methods  which  will  be  discussed. 

It  is  possible  using  a reverse  flow  theorem  to  determine  the  force  which 
is  normal  to  the  fin  planform  and  acts  on  the  fin-body  conbination.  We 
then  multiply  the  fin-body  normal  force  by  the  ratio,  (K^  + Kg),  to 
obtain  an  estimate  of  the  normal  force  acting  on  the  fin.  We  divide  the 
normal  force  acting  on  the  fin  by  the  fin-alone  normal  force  for  unit 
angle  of  attack  and  obtain  an  equivalent  angle  of  attack  (^“eq)v,i 
the  i'th  fin  due  to  the  vortex.  This  procedure  can  be  carried  out  using 
either  linear  theory  or  slender-body  theory  provided  the  same  theory  is 
used  throughout.  Any  bias  is  minimized  since  a normal-force  ratio  is 
utilized  to  obtain  an  equivalent  angle  of  attack. 

The  next  step  is  to  use  the  nonlinear  normal-force  curve  for  the 
fin  alone  to  determine  the  fin  force  associated  with  (^%q^,i* 
fin  force  depends  at  what  point  on  the  normal-force  curve  the  fin  is 
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operating  before  account  is  taken  of  (Aa^q)^  The  associated  force 
acting  on  the  body  is  then  obtained  by  multiplying  the  fin  normal  force 

by  the  ratio,  K3/K^^.  The  procedure  can  be  applied  to  all  four  fins  to 

obtain  the  normal-force  and  side-force  components  of  the  fin-body 
combination. 

With  respect  to  rolling  moment,  it  is  possible  to  obtain  the 
rolling  moment  acting  on  the  wing-body  due  to  the  fin  twist  associated 
with  a vortex  using  reverse  flow  theorems.  This  rolling  moment  is  due 
in  part  to  direct  roll  associated  with  the  fin  in  question  and  by  reverse 
roll  associated  with  the  other  fins.  If  the  rolling  moment  (about  the 
body  axis)  is  divided  by  the  fin-alone  normal  force  calculated  by  the 
same  theory  as  the  rolling  moment,  then  an  effective  lever  arm  for  the 
fin  is  established.  This  lever  arm  includes  the  effects  of  reverse  roll 

and  is  thus  closer  to  the  body  axis  than  if  reverse  roll  had  been 

neglected.  The  nonlinear  rolling  moment  due  to  the  fin  can  then  be 
calculated  by  multiplying  the  nonlinear  fin  normal  force  by  this  lever 
arm. 

The  detailed  derivations  of  the  equations  based  on  this  method  are 
now  given. 

B.3.  BOUNDARy  CONDITIONS  INDUCED 
ON  FINS  BY  VORTICES 

The  procedure  for  calculating  the  distribution  of  induced  angle  of 
attack  normal  to  the  fins  is  as  follows: 

1.  Specify  the  spanwise  fin  coordinates  at  which  the  induced 
angle  of  attack  is  to  be  calculated. 

2.  Calculate  the  velocities  induced  at  those  coordinates  by  the 
Biot-Savart  law  for  all  the  vortices  present. 

3.  From  the  velocities  computed  for  each  point,  determine  the 
component  normal  to  the  fin. 

4.  Depending  upon  the  option  of  the  user,  limit  the  velocity 
which  may  be  induced  on  the  fin. 

The  code  gives  the  user  a nunber  of  options  concerning  the  induced  angle 
of  attack.  The  options  are:  (1)  impose  no  limit  on  the  angle  of  attack 

which  may  be  induced  on  the  fin;  (2)  assume  that  each  vortex  has  a 
viscous  core;  or  (3)  limit  the  maximum  velocity  which  may  be  induced  on 
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a fin  to  a specified  fraction  of  V^.  If  either  option  (1)  or  (2)  is 

selected,  a is  set  equal  to  the  arc  tangent  of  the  ratio  of  the 

^ • -1 
induced  velocity  to  V^.  If  option  (3)  is  selected,  ^ xs  set  equal 

to  the  ratio  of  the  induced  velocity  to  or  whichever  is 

smaller. 

B.4.  DETERMINATION  OF  LOADING 
ACTING  ON  FINNED  SECTION 

Since  xt  is  based  on  integral  theorems,  the  reverse  flow  method 
cannot  give  the  loading  on  individual  fins.  Rather,  it  gives  the  overall 
loading  on  the  finned  section  due  to  the  flow  field  seen  by  each  fin 
This  is  sufficient  for  small  angles  of  attack  for  which  the  fin  loading 
characteristics  are  linear.  However,  when  the  fin  loading  characteristics 
are  nonlinear,  an  extension  to  the  reverse  flow  method  is  required.  The 
extension  used  in  the  present  method  is  based  on  the  equivalent  angle  of 
attack  concept.  The  procedure  we  will  use  is  as  follows; 

1.  Compute  the  normal  force  or  side  force  acting  on  the  fin-body 
configuration  due  to  the  fin  twist  associated  with  the  vortices  using  a 
combination  of  reverse  flow  theory  and  slender-body  theory. 

2.  Estimate  the  fraction  of  the  normal  force  or  side  force  acting 
on  the  fin  by  multiplying  the  quantity  calculated  in  (1)  by  the  ratio, 

V (%  + %)  • 

3.  Divide  the  fin  force  of  (2)  by  the  slope  of  the  slender-body 
theory  fin-alone  normal-force  curve  to  obtain  an  equivalent  angle  of 
attack,  (Aa„„)„.  Using  this  quantity  with  the  fin-alone  nonlinear  normal- 
force  curve  will  give  the  nonlinear  fin  normal  force  in  the  presence  of 
the  body. 

4.  Using  the  same  theory  as  in  (2)  above,  calculate  the  total 
rolling  moment  acting  on  the  finned  section  due  to  the  fin  twist 


associated  with  the  vortex. 

5.  Divide  the  result  of  (4)  by  the  result  of  (2)  to  get  a moment 
arm  for  the  force  due  to  (Aa^g)^. 
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B.4.1  Computation  of  Equivalent  Angle 
of  Attack  for  Each  Fin 


Consider  now  a cruciform  finned  section  with  a circular  body  in 
direct  flow  and  in  reverse  flow  as  shown.  We  will  consider  the  force 
normal  to  the  plane  of  fins  2 and  4 first  and  generalize  the  results 
to  an  arbitrary  fin.  The  fins  are  undeflected. 


Fin  4 Fin  4 


Sketch  B.l,  Direct  and  reverse  flow. 


Let  Sg  be  the  body  planform  area  and  let  Sj^,  , s^,  and  be  the 
planform  areas  for  the  fins.  Let  and  be  the  fin  loading  coeffi- 

cient and  the  local  angle  of  attack  in  direct  flow,  and  let  P and  a 
be  the  corresponding  quantities  in  reverse  flow.  Then,  the  reverse  flow 
theorem  (ref.  9,  p.  221)  gives 


/ / ‘■r 

Sg+S^+S^ 

+S3,S^ 


Sg+S^+S 


2 


+s„+s 

3 


4 


P ,a  dS 
d r 


(B.l) 


Choose  the  following  specific  values  of  angles  of  attack  and  pressure 
coefficients 


Direct 

Reverse 

“d 

0 on 

a 

r 

= 1 

on 

II 

0° 

a . on 
v,i 

s^.s. 

a 

r 

= 1 

on 

p.  = 

P. 

a 

= 0 

on 

S,  ,s^ 

Ihen,  (Brl)  becomes 


/ / VS  " /[%,aVS  + //  %,4^3dS  (B.: 

Sb+Si-^S^  S4 

where  P3  is  the  loading  for  a wing-body  combination  at  unit  angle  of 
attack  and  no  yaw  or  pitch  control.  The  force  normal  to  the  plane  of 
fins  2 and  4 acting  on  the  finned  section  is 


r 


T 


SJ^S^+S 


PddS 


(B.3) 


2 4 

Substituting  (B.2)  into  (B.3)  gives 


s„  t.e. 
• m 


s„  t.e. 
m 


a i.e.  a I.e. 


where  the  dummy  variable  t refers  to  integration  along  the  span  of  the 
fins.  Let  the  span-load  distribution  associated  with  P^  be  (cc^)g 
such  that 


(cc^). 


t.e. 


:.e. 


PgdX 


(B.5) 


Now  we  assume  that  ^ varies  slowly  in  the  chordwise  x-direction  so 
that  substitution  of  (B.5)  into  (B.4)  gives 

r^m 

N = ‘^oo  J J <S.€ 

a a 

The  normal  forces  acting  on  Sg  and  are  not  given  by  this  result 
since  part  of  the  load  due  to  each  integral  can  appear  on  the  body  and 
the  opposing  fin.  Here  (cc^)^  is  the  span  loading  on  the  fins  for  the 
configuration  in  reverse  flow  with  both  body  and  panels  at  unit  angle 
of  attack  with  no  yaw  or  pitch  control.  From  reference  9,  equation 
(5-15),  we  have  the  slender-body  theory  expression  for  the  span  loading: 
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(CC_^)3 


(B.7) 


In  order  to  obtain  the  normal  force  on  each  individual  fin  and 
thereby  carry  out  item  2 of  the  procedure  described  above,  we  multiply 
equation  (B.6)  by  K^/(K^  + K^)  and  consider  each  term  to  be  the  force 
acting  on  the  fin  over  which  the  integration  is  taking  place.  We  can 
follow  the  same  procedure  for  fins  1 and  3 with  the  following  result  for 
fin  i: 

*^fin  i = +■  Kg  J “v,i<^)  (cc^>3dt  (B.8) 


p 

The  slender-body  normal-force  slope  for  the  frn  is  irq^(Sj^-a)  . Hence, 
carrying  out  item  3 of  the  procedure  described  at  the  beginning  of  this 
section,  we  find 


(Aa  ) = 

eq  V 


fin  i 


cxv,i(t)  (cc^)3dt 


(B.9) 


Without  the  vortex  present,  the  i'th  fin  operates  at  (a  ^ • The 

difference  in  normal  force  at  p + ^°^eq  p^  i 

the  fin-alone  nonlinear  normal-force  curve  is  the  nonlinear  normal  force 
on  the  fin  associated  with  the  vortex.  It  must  be  multiplied  by  (Kg/K^) 
to  obtain  the  nonlinear  contribution  due  to  body  carryover. 


B.4.2  Computation  of  the  Moment  Arm  for 
the  Rolling  Moment  Associated  with 
Each  Fin  due  to  the  Vortex 

The  reverse-flow  theorem  for  all  the  fins  with  the  body  aligned  in 
the  free-strecim  direction  (i.e.,  a = 0 on  body)  is 


+S3+S 


4 


adPrdS 


afPddS 

Si+S^ 

+S3+S^ 


(B.IO) 
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where  t is  the  spanwise  coordinate  for  each  fin  measured  from  the  body 
axis.  Note  that  corresponds  to  the  distribution  of  angle  of  attack 

which  would  be  seen  by  a configuration  rolling  at  the  rate  of  cjo  radians 
per  second  clockwise.  Then,  (B.IO)  becomes 


has  the  sign  convention  shown  in  the  sketch  below 


Sketch  B.2.  Sign  convention  for  a, 


The  total  rolling  moment  for  the  finned  section  is 


Substituting  (B.ll)  into  (B.12)  gives 


M 


I CJD 


s_  t.e. 
• m 


s t.e. 

-m  - 


- a i.e.  a J.e. 


P„dx  dt 

t? 


/ 


m 


t.e. 

• m 


a i.e. 


j j «V,4(^.t) 

a 2.e. 


P„dx  dt 


(B.13) 


Let  the  span  loading  for  the  reverse  flow  case  above  be 

t.e. 


V 


(cc.)  = — 

re  CO 


P dx 


(B.14) 


.e. 


Recall  our  original  assiomption  that  a . is  constant  in  the  chordwise 
direction,  and  substitute  (B.14)  into  (B.13).  The  result  is 


^ 


p - 
- / a 

J 


(t)  (cc^)^dt  + 


^(t)  (cc_g)^dt 


L a 


/m 

- J %,4 


m 


(t) (cc J dt 

Zf  6 


(B.15) 


The  integrals  in  equation  (B.15)  do  not  represent  the  actual  rolling 
moments  exerted  on  the  individual  fins.  Rather,  they  represent  the  con- 
tribution to  the  total  rolling  moment  due  to  the  distribution  of  angle 
of  attack  induced  on  each  fin.  To  determine  the  actual  rolling  moment 
for  each  fin,  we  must  determine  the  effect  of  each  fin  on  every  other 
fin.  Let  us  define  the  following  quantities. 


M 


s rolling  moment  induced  on  fin  i 


(B.16) 


fin  i 


r - 

I \r  1 


(t)  (cc  J dt 


(B.17) 


-1  , i - 1 
+1  , i - 2 
+1  , i - 3 
-1  , i - 4 


= fraction  of  L contributing  to  rolling 
^ moment  of  fin  i 


(B.18) 


(B.19) 


Using  the  above  definitions,  we  can  write 


Note  that 


11  12  13  14 


21  22  23  24 


31  32  '^33  ^^34 


41  42  43  44 


k . + k . + k . + . = 1 

11  21  31  41 


(B.20) 


so  that  equation  (B.15)  is  recovered. 

In  principle,  the  k^j  can  be  computed  by  slender-body  theory  or 
by  linear  theory  computer  codes  such  as  CRFWBD  (Appendix  F) . However, 
for  i ^ j,  the  coefficients  are  generally  considerably  less  than  one 
although  not  negligible.  However,  we  neglect  them  in  the  present  approxi- 
mate method  for  calculating  the  individual  fin  rolling  moments  and  assume 


r 1 , i “ j 
“lo  . i / j 


(B.21) 


Substituting  the  approximation  of  (B.21)  into  (B.20)  gives  the  following 
rolling  moment  for  each  fin; 


\in  i “ 
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The  slender-body  theory  expression  for  (cc.)  is  given  by  Adcuns 

Jc  & 

and  Dugan  (ref.  B.2,  eq.  (18)): 


cos  29  tanh~^  ( ^ tanh  ^ 

TT  \ sin  2y  / 


( tan  29  \1 
\tan  2T/J 


^^[K(k^)sin  40  - 2k^  cos  A^K (k^ ) Z (A^  ,k^  )] 


where 


(B.23) 


(B.24) 


t*' 


cos  2Y  = — — 

2R® 

k^  = sin  27 

tt/2 

dz 

1 - k^  sin^  z 
o 1 


K(k^)  = r 


(B.25) 


i 

(B.26) 

I 

(B.27)  I 


= complete  elliptic  integral 
of  first  kind  with 

modulus  k^  (B.28) 


A 

1 


(B.29) 


Z(A^,k^)  = Jacobi  zeta  function 


(B.30) 


K(k^)Z(A^,kJ 


K(k^)E(A^,k^)  - E(Tr/2,kjF(A^,k^) 


(B.31) 


F(A^,kJ 


dz 

'\l  1 - k^  s in-  £ 


incomplete  elliptic  integral 
of  the  first  kind 


(B.32) 
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E(7T/2,k^) 


1 ~ sin^  z dz  = complete  elliptic  integral  (B.33) 
^ of  second  kind 


E(A^,k^) 


1 - k^  sin^  z dz  = incomplete  elliptic  integral  (B.34) 
^ of  second  kind 


In  equation  (B.23)  as  6 — y;  that  is,  as  t — a,  the  first  of  the  two 
principal  terms  becomes  indeterminate.  It  is  possible  to  show  that 


Lim  cos  20  tanh"^  (sin~2r)  " tanh"^  ^ 


tan  20  \ 

tan  2y  Jj 

-cos  2Y  log(cos  2Y)  (B.35) 


Thus,  the  span  loading  at  the  juncture  between  the  fin  and  body,  t = a. 


(cc.)  = r -cos  2Y  log(cos  2Y)  + 0.5  K(k,  )sin  40] 

i 6 TT  L 1 -I 


(B.36) 


To  obtain  the  moment  arm  for  the  fin  load  represented  by  (Aa^q)^, 


we  divide  M 


e i 5 that  is, 

^fin  i ^ 


f 

^fin  i 


(B.37) 


The  total  rolling  moment  is  then  obtained  by  multiplying  Yv^  by  each 
corresponding  fin  nonlinear  normal  force  as  obtained  in  item  3 of  the 
procedure  described  above. 
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APPENDIX  C 


PROCEDURE  FOR  COMPUTING  STRENGTHS  AND 
POSITIONS  OF  FIN  TRAILING  VORTICES 


in  setting  up  a model  of  a single  vortex  per  fin  we  must  specify 
the  vortex  strengths  and  the  lateral  position  of  the  center  of  vorticity. 
Using  the  equivalent  angle  of  attack  concept,  we  can  obtain  a good  esti- 
mation of  the  fin  normal  force  and  lateral  center  of  pressure  position. 
The  single-vortex-per-panel  model  is  compatible  with  this  input  data. 

To  go  to  more  vortices  per  fin  requires  knowledge  of  the  span-load 
distribution. 

At  low  angles  of  attack,  the  span-load  distribution  of  slender  wings 
tends  to  be  elliptical.  At  high  angles  of  attack,  the  wing  loading  tends 
to  be  uniform,  and  for  a delta  wing  the  span  loading  tends  to  be  linear. 
As  an  approximation  to  actual  span  loadings  consider  a superposition  of 
an  elliptical  loading  plus  a linear  loading.  We  will  now  derive  expres- 
sions for  the  normal  force  and  rolling  moment  associated  with  elliptic 
and  linear  span  loadings. 


T 

m 


(C.l) 


By  definition. 


(cc_^)e  = a^Vl  - 


s 

I 


m 


(cc  Jdy 


CN  = 


rm 

/ y(cc^)dy 


CRM  = 


^R^r 


(C.2) 


(C.3) 


(C.4) 
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Now  the  first  integral  can  be  written  as 


J (cc^)dy  = {s^  - a)  J (cc^dT] 


= (s  - a) 
m 


X J. 

a^  J a - T])dT)  + J Vl  - 


rf  dT) 


(C.5) 


f 


to 


J = (s^  - a) 


(C.6) 


Thus 


CN  = 


3 - a 

m 


For  the  rolling  moment,  we  obtain 


/m  f-m 

(cc^)dy  + 


CRM  = 


®R^r 


a(Sm  - a)  /a. 


(C.7) 


®R^r 


X 1 

J T](l  - T])dT)  + \ - ri^  dTi 


a(Sm  - a)  / a. 


Sr^ 


®R^r 


(C.  8) 


The  lateral  location  of  the  center  of  pressure  is  simply  given  by 

t.cm  (^^-x) 


CN 


a + (s  - a) 
m ' 


(x  • I *e) 


(C.9) 
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or  in  nondimens ional  fom 


■Y  = 
s - a 
m 


- 1 s) 


(C.IO) 


Solving  for  find 


cly  1 TT 

■>->  T ” /I 


1 - 1 
2 ^ ~ 6 
1 TT  - 
3 ■ 4 


(C.ll) 


This  gives  the  ratio  ^e/^L  terms  of  the  lateral  location  of  the 
center  of  pressure.  The  fin  normal-force  coefficient,  CN,  then  permits 
us  to  obtain  aj^  frora  eg’jis.tion  (C.  7)  as  follows: 


S^CN 


(s  - a) 
' m ' 


(^5^) 


(C.  12) 


Having  a^  and  corresponding  to  the  data,  we  now  must  find  the 
vortex  strength,  r,  and  its  lateral  position,  y*.  The  lateral  position, 
y*,  is  at  the  centroid  of  the  trailing  vorticity  given  by  the  following 
expression  (see  ref.  9^  pp.  144-153) : 


I 


m d(cc.) 

y dy 


® 3 + 


/ 


'm  d (cc^) 
dy 


dy 


(C.13) 


Carrying  out  the  indicated  operations  yields 


y*  - a 
•'^v  

s - a 
m 


, + 7 a 
2 4 e 

a,  + a 


(C.14) 


Assuming  that  the  part  of  the  bound  vortex  outside  the  body  represents 
the  fin  normal  force,  we  can  write 


p V r (y*  - a)  “■ipv^S  C 
00  'Jfy  2 ^00  00  '“N 


(C.15) 
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Substituting  (C. 14)  into  (C.16)  gives 


APPENDIX  D 


MEIHOD  FOR  COMPUTING  VORTEX  POSITIONS 
AND  STRENGTHS  OVER  AFTERBODY  SECTION 

D.l.  INTRODUCTICai 

The  purpose  of  this  appendix  is  to  describe  the  method  used  by  the 
code  MLQADS  to  track  vortices  downstream  of  the  first  finned  section. 
The  method  consists  of  three  procedures:  (1)  determination  of  the 

starting  positions  of  the  afterbody  vortices;  (2)  tracking  of  all  the 
vortices  present  over  the  afterbody  section;  and  (3)  determination  of 
the  strengths  of  the  two  afterbody  vortices  as  functions  of  axial 
location.  The  code  MLOADS  uses  the  method  described  in  reference  3 for 
procedure  (1)  with  one  exception.  In  the  method  of  reference  3,  after- 
body vortices  are  not  started  unless  the  direction  of  the  crossflow 
velocity  vector  is  within  90°  of  the  vertical  eixis,  z^.  In  the  present 
version  of  the  code,  the  user  selects  an  angle  OMEGA.  If  the  crossflow 
velocity  vector  is  not  within  +0MEGA  of  z^,  afterbody  vortices  are  not 
allowed  to  form.  Procedures  (2)  and  (3)  are  described  in  the  sections 
below. 

D.2.  TRACKING  OF  VORTICES  OVER 
THE  AFTERBODY  SECTION 


We  assume  that  the  vortex  centers  move  as  Lagrangian  fluid  points 
in  a crossflow  plane  moving  with  the  speed  V^cos  a^.  The  complex 
velocity  for  vortex  j in  the  presence  of  a circular  body  and  the 
other  vortices  may  be  obtained  from  equation  (A. 10)  by  setting 
dv/da  = 1 and  r^  = a.  The  result  is 

n 

i E '’k 

k=l 

Mj 


+ 


i 

2r 


Oj  - a»/3j 


(D.U 


where  n is  the  number  of  vortices  present.  We  may  now  repeat  the 
development  following  equation  (A. 12)  to  get 
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and 


(D.2) 


1 

cos  a 

00  Q 


1 

cos  a 

00  Q 


Again,  the  present  version  of  the  code  uses  small  angle  approximations 
and  cos  in  (D.2)  is  set  equal  to  one.  The  differential  equations 

(D.2)  and  the  algebraic  relations  (D.l)  provide  a system  of  ordinary 
first-order  initial  value  differential  equations  which  can  be  solved  by 
standard  integration  routines.  The  system  is  complete  if  no  afterbody 
vortices  are  present.  If  afterbody  vortices  are  present,  two  of  the 
are  unknown  functions  of  x.  Differential  equations  for  the  afterbody 
vortex  strengths  are  derived  in  the  next  section. 


D.3.  DERIVATION  OF  DIFFERENTIAL  EQUATIONS 
FOR  THE  STRENGTHS  OF  THE 
AFTERBODY  VORTICES 


Consider  a crossflow  plane,  shown  in  sketch  D.l,  containing  the 
afterbody  cross  section,  and  let  afterbody  vortices  F^  and  F„  appear 

Li  H 

at  stations  a and  a in  the  crossflow  plane.  The  other  vortices  are 

Xj  h 

not  shown  for  simplicity. 


Sketch  D.l.  Computational  plane. 


In  a small  eixial  distance,  dx,  the  afterbody  increments  in  normal  force 
and  side  force  are  given  by  the  vortex  impulse  theorem  (eq.  (28)). 


dN  - i dY  » p V^d 
' 00  00 


(t),3) 
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If  r,  and  r„  are  constant,  then  the  resulting  terms  on  the  right-hand 

li  R 

side  of  equation  (D.2)  are  differentials  of  vortex  position  and  are 


associated  with  inviscid  forces  on  the  afterbody.  If  r,  and  r„  are 

Xj  I\ 


not  constant  because  of  vortex  feeding  sheets,  then  additional  terms 


containing  differentials  of  and  r„  appear  which  are  associated  with 

Xj  R 


crossflow  drag.  If  we  use  the  subscript  v to  indicate  that  part  of 
dN  - i dY  due  to  viscous  crossflow,  then 


- i <«v-  f.''.  [ ("i,  - + (“r  - ^)  ■’i-r] 


(D.4) 


NOW  we  will  assume  that  the  crossflow  drag  is  in  the  seune  direction 
as  the  crossflow  velocity  vector,  9.  Then, 


dN  - i dY  = D (sin  0 - i cos  0) 

V V V ' 


(D.5) 


where  is  the  crossflow  drag  force  on  the  elemental  length  of 


cylinder.  Ax.  In  terms  of  the  crossflow  drag  coefficient,  is 


= c , q^sin^a  ^ 2a  dx 
V d ^ cr 

c 


(D.6) 


Siibstitution  of  (D.5)  and  (D.6)  into  (D.4)  gives 
PcovJ  (a,  - + (a„  - 2a  c,  q^  sin 


\( 

) dx 

■ i:- 

) dx 

(D.7) 


Equation  (D.7)  gives  two  scalar  first-order  initial  value  ordinary 
differential  equations  which  can  be  solved  simultaneously  with 
equation  (D.2). 
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APPENDIX  E 


MLOADS,  A CCMPUTER  PROGRAM  FOR  CALCULATING 
THE  COMPONENT  FORCES  AND  MOMENTS  OF 
ARBITRARILY  BANKED  CRUCIFORM  MISSILES 
WITH  CONTROL  DEFLECTIONS 


E.l.  INTRODUCTION 


The  purpose  of  this  appendix  is  to  describe  the  operation  of  the 
computer  code  MLOADS  in  sufficient  detail  to  permit  understanding  and 
use  of  the  program.  The  program  computes  forces  and  moments  for  each 
fin,  for  each  section  of  the  configuration,  and  for  the  complete  config- 
uration. Minimum  drag  is  not  computed.  Intermediate  output  gives  vortex 
strengths  and  positions  along  the  body.  The  code,  through  the  integer 
control  variable, NFIN,  is  capable  of  handling  the  following  missile 
configurations; 


NFIN 

NFIN 

NFIN 

NFIN 


1:  body-alone  consisting  of  nose  and  cylindrical  afterbody 

2:  body-tail  combination  (f inner) 

3;  body-canard  combination  (body-wing) 

4;  body-canard-tail  conibination  (body-wing-tail) 


Typical  exeunples  of  the  above  configurations  are  given  in  the  sketch 
below. 


PrtctdKg  pad  Hank 
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Body-canard-tail 


Sketch  E.l.  Excimples  of  missile  configurations. 

A description  of  the  engineering  method  used  is  given  in  reference  3 and 
in  the  main  text  and  Appendices  A through  D of  this  report. 

The  program  is  written  in  FORTRAN  IV  (029  punch)  and  has  been  run 
on  the  CDC  6600  and  7600  machines.  No  files  other  than  INPUT  and  OUTPUT 
are  required.  The  program  conforms  to  the  ANSI  standard  Fortran  language 
except  for  the  program  definition  statement  and  for  several  DATA  state- 
ments. For  a body-tail  combination,  the  running  time  is  roughly  0.08  CPU 
second  per  data  point  on  the  CDC  7600.  For  a body-canard-tail  combination 
the  ruj'ning  time  is  about  0.5  CPU  second  per  data  point  on  the  CDC  7600. 
Run  times  on  the  CDC  6600  are  about  six  times  those  for  the  CDC  7600. 

E.2.  PROGRAM  DESCRIPTION 

For  computational  purposes,  the  missile  to  be  studied  is  divided 
into  four  sections  as  shown  in  figure  1 of  the  main  text.  The  nose  sec- 
tion is  defined  to  be  from  the  nose  tip  to  the  leading  edge  of  the  root 
chord  of  the  first  finned  section.  The  canard  section  (first  finned 
section)  is  defined  to  be  from  the  end  of  the  nose  section  to  the  trailing 
edge  of  the  first  set  of  fins.  The  afterbody  section  is  defined  to  be 
from  the  end  of  the  canard  section  to  the  leading  edge  of  the  root  chord 
of  the  second  finned  section.  The  tail  section  (second  finned  section) 
is  defined  to  be  from  the  end  of  the  afterbody  section  to  the  trailing 
edge  of  the  second  set  of  fins.  The  computations  proceed  as  shown  in 
figure  E.l.  The  branching  for  the  four  possible  configurations  is  shown 
in  figure  E.2.  The  main  program  MLOADS  controls  the  flow  of  operations. 

It  calls  the  routines  which  control  the  operations  for  each  section  of 
the  missile,  A more  complete  description  of  the  program  operation  is 
given  below. 


INPT 

1.  X,  y for  Canard  Fins 

2 . , X,  y for  Tail  Fins 

3 . Linear  Interference  Parameters  (Optional 
4 . Body  Geometry 

5 . Flight  Conditions 

6 . Run  Control  Parameters 

t -1 

NOSE 

Computes 

1.  Nose  Loads 

2.  Positions  of  Nose  Vortices  Using 

Empirical  Correlations 

3.  Strengths  of  Nose  Vortices  Using  c^j 

and  Vortex  Impulse  Theorem 

i 1 

CANARD 

Computes 

1.  Positions  of  Nose  Vortices  over 

Canard  Section 

2.  Individual  Fin  Loads  by  Equivalent 

Angle  of  Attack  Method 

3.  Strengths  and  Positions  of  Canard 

Trail inq  vortices 

i 

AFTBOD 

Computes 

1.  Positions  of  Vortices  over  Afterbody 

as  a Function  of  Axial  Position 

2.  Strengths  of  Afterbody  Vortices 

3.  Loads  on  Afterbody  Section 

♦ 

TAIL 

Computes 

1.  Positions  of  Vortices  over  Tail 

Section 

2.  Individual  Fin  Loads  by  Equivalent 

Angle  of  Attack  Method 

PRINT  INPUT  DATA 


PRINT  NOSE  L07U5S 
‘AND  VORTEX  VALUES 


PRINT  INDIVIDUAL  FIN 
LOADS,  TOTAL  LOADS 
-FOR  CANARD  SECTION, 
AND  VORTEX  STRENGTHS 
AND  POSITIONS 


PRINT  POSITIONS  AND 
STRENGTHS  OF  AFTER- 
BODY VORTICES  AS  A 
FUNCTION  OF  AXIAL 
POSITION  AND  TOTAL 
LOADS  ON  AFTERBODY 
SECTION 


PRINT  INDIVIDUAL  FIN 
'LOADS  AND  TOTAL  LOADS 
FOR  TAIL  SECTION 


PRINT  SUMMARY  OF  TOTAL  MISSILE  FORCES  AND  MOMENTS 


Figure  E.l.  Flow  of  progreim  computations. 
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E.2.1  Calculation  Procedure 


After  reading  in  the  run  identification  information,  the  subprogrcun 
INPT  reads  in  body  geometry,  run  control  parameters,  flow  conditions,  emd 
any  required  empirical  information  not  built  into  the  program.  All 
length  quantities  are  then  normalized  by  the  body  radius  and  all  areas 
are  normalized  by  the  square  of  the  body  radius.  The  user  has  the  option 
of  allowing  the  code  to  compute  the  necessary  interference  pareuneters  or 
of  reading  them  in.  If  the  user  selects  the  former  option,  the  methods 
of  slender-body  theory  and  reference  10  are  employed.  The  user  may  read 
in  any  or  all  of  the  interference  parameters. 

The  user  must  read  in  the  wing-alone  characteristics  for  each  set 
of  fins.  These  characteristics  include  the  normal-force  coefficient, 
spanwise  location  of  the  fin  center  of  pressure  and  chordwise  location 
of  the  fin  center  of  pressure  as  functions  of  angle  of  attack. 

After  all  the  input  has  been  read  in  and  all  required  interference 
factors  have  been  computed,  the  main  progreun  MLOADS  calls  subroutine  NOSE. 
This  routine  first  computes  the  axial  starting  location  of  the  nose  vor- 
tices (if  present) . The  locations  and  strengths  of  the  nose  vortices  in 
the  crossflow  plane  at  the  leading  edge  of  the  root  chord  of  the  first 
finned  (canard)  section  are  computed  next.  The  siibroutine  then  computes 
the  forces  and  moments  on  the  nose  and  returns  control  to  MLOADS. 

If  NFIN  = 1,  the  driver  program  prints  a summary  of  the  overall 
forces  and  moments  and  reads  in  a new  data  set.  If  NFIN  ^ 2,  subroutine 
CANARD  is  called.  This  routine  computes  the  fin  and  body  loadings  for  the 
first  set  of  fins  downstream  of  the  nose.*  The  individual  fin  forces 
and  moments  in  the  absence  of  nose  vortices  are  computed  first.  The 
routine  then  calls  the  routine  CRUTRJ  which  computes  the  positions  of 
the  nose  vortices  over  the  canard  section.  The  discrete  vortex  model 
used  to  represent  the  nose  vorticity  may  not  adequately  represent  the 
flow  when  the  separated  flow  from  the  nose  is  diffuse  or  when  a vortex 
is  close  to  a fin.  With  the  present  model,  it  is  possible  for  a vortex 
to  be  captured  by  its  image  and  to  travel  completely  around  the  body. 

Such  behavior  has  not  been  observed  in  flow  visualization  studies  and 
the  predicted  loads  obtained  when  a vortex  is  "captured"  do  not  agree 

*Note  that  if  the  configuration  is  a body-tail  combination,  the  tail 
section  results  are  computed  by  subroutine  CANARD  and  are  output  under 
the  headings  for  "CANARD  SECTION." 


with  experiment.  We  have  found  that  the  present  model  gives  best  results 
when  the  presence  of  the  firs  is  ignored  during  the  vortex  tracking 
procedure.  The  user  can  nxercise  this  option  by  setting  the  logical 
variable  NEARC  equal  to  .TRUE.  This  option  is  recommended.  The  user 
may  also  hoose  to  ignore  the  nose  vortices  over  all  or  part  of  the 
missile  if  he  believes  they  have  dispersed.  This  control  is  achieved 
through  the  parameter  NVORT. 

When  CRUTRJ  has  finished  tracking  the  nose  vortices,  it  stores  their 
locations  in  the  crossflow  planes  at  the  fin  area  centroid  and  at  the 
fin  trailing  edge  and  returns  control  to  CANARD.  Subroutine  CANARD  then 
calli.  REVFLO  which  computes  the  equivalent  angle  of  attack  and  spanwise 
location  of  the  center  of  pressure  for  the  loading  due  to  the  nose  vor- 
tices. Subroutine  REVFLO  assumes  the  vortices  to  be  infinite  line  vor- 
tices parallel  to  the  body  axis.  The  vortex  locations  at  the  fin  area 
centroid  are  used. 

When  REVFLO  returns  control  to  CANARD,  the  total  fin  forces  and 
moments  and  body  forces  and  moments  (for  the  canard  section)  are  computed. 
The  strengths  and  positions  of  the  vortices  shed  from  each  fin  are  calcu- 
lated next.  After  the  above  information  has  been  printed  out,  control  is 
returned  to  MLOADS. 

If  NFIN  = 2,  the  driver  program  prints  a summary  of  the  overall 
forces  and  moments  and  reads  in  a new  data  set.  If  NFIN  > 3,  sub- 
routine AFTBOD  is  called.  This  routine  tracks  the  trailing  vortices 
from  the  first  set  of  fins  and  the  nose  vortices  (if  present)  over  the 
afterbody  section.  The  routine  also  decides  if  afterbody  vortices  should 
be  shed  and  computes  their  starting  locations.  If  afterbody  vortices  are 
formed,  they  are  tracked  along  with  the  other  vortices  present.  The 
strengths  of  the  afterbody  vortices  as  functions  of  axial  position  are 
also  computed.  Once  the  vortex  trajectories  have  been  computed  for  the 
afterbody  section,  the  forces  and  moments  on  the  afterbody  section  are 
computed  and  control  is  returned  to  MLOADS. 

If  NFIN  = 3,  MLOADS  prints  a summary  of  the  overall  forces  and 
moments  and  reads  in  a new  data  set.  If  NFIN  = 4,  subroutine  TAIL  is 
called.  This  routine  proceeds  in  a manner  similar  to  CANARD  except  that 
the  trailing  vortices  of  the  tail  fins  are  ignored.  When  TAIL  has 
finished,  it  returns  control  to  MLOADS.  Subroutine  MLOADS  then  computes 
and  prints  out  the  total  forces  and  moments  for  the  configuration  and 
reads  in  a new  data  set. 
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£•2.2  Program  Limitations  and  Precautions 


i 


The  program  makes  a number  of  assumptions  about  the  missile  configu- 
ration  and  the  flow  field  which  are  listed  in  section  6 of  the  main  text.  i 

i 

E.2.3.  Description  of  Subroutines  ; 

Subroutine  INPT  reads  in  all  input  and  nondimensionalizes  it. 

Subroutine  NOSE  computes  nose  vortex  positions  and  strengths  and 
nose  forces  and  moments. 

Subroutine  CANARD  controls  the  subroutines  which  compute  vortex 
strengths  and  positions,  individual  fin  forces  and  moments,  body  forces 

and  moments  and  total  forces  and  moments  for  the  canard  section.  ' 1 

Subroutine  ALFAQ  computes  equivalent  angles  of  attack  for  individual 
fins  in  the  absence  of  vortices.  j 

Siabroutine  VORTEX  computes  the  strengths  and  positions  of  the  vortices  1 

shed  from  the  trailing  edges  of  the  canard  fins. 

Sxibroutine  CRUTRJ  computes  the  paths  of  free  vortices  in  the  presence 
of  a cruciform  wing-cylindrical  body  section. 

Subroutine  FCT  computes  derivatives  needed  by  the  differential 
equation  solver  of  CRUTRJ. 

Subroutine  SHAPE  provides  the  body  and  fin  coordinates  needed  by  FCT. 

Subroutine  DASCRU  is  the  differential  equation  solver  used  by  CRUTRJ 
and  AFTBOD. 

Subroutine  REVFLO  computes  the  equivalent  angles  of  attack  and  span- 
wise  locations  of  the  centers  of  pressure  for  a set  of  fins  due  to  the 
presence  of  vortices. 

Subroutine  SIMSON  is  a Simpson’s  rule  integration  package  used  by 
REVFLO. 

Subroutine  VEL  computes  the  velocities  induced  on  a set  of  fins  by 
the  vortices  present  for  use  by  REVFLO. 

i Subroutines  CELl,  CEL2,  ELIl,  and  ELI2  are  IBM  SSP  subroutines  for 

1 computing  elliptic  integrals.  They  are  called  by  REVFLO. 

I ; Subroutine  LNTRP  is  a linear  interpolation  routine. 
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Subroutine  CURVES  is  a BLOCK  DATA  routine  for  initializing  a number 
of  empirical  tables  used  in  the  progreim. 

Subroutine  AFTBOD  computes  vortex  positions  and  strengths  and  body 
loading  for  the  afterbody  section. 

Subroutine  DERIV  computes  derivatives  needed  by  the  differential 
equation  solver  of  AFTBOD. 

Svibroutine  TAIL  controls  the  subroutines  which  compute  vortex 
strengths  and  positions,  individual  fin  forces  and  moments,  body  forces 
and  moments  and  total  forces  and  moments  for  the  tail  section. 

Subroutines  CHRT8,  EQ24,  EQ26,  EQ30,  EQ31,  and  CH1416  are  used  to 
compute  Kg  and  Xg  based  on  the  methods  of  reference  10.  These 
subroutines  are  described  in  reference  E.l. 

Subroutine  INTFAC  calculates  any  interference  parameters  which  have 
not  been  input. 

E.3.  DESCRIPTION  OF  INPUT 


This  section  describes  the  input  for  the  program  MLOADS.  In  the 
following  discussion,  the  contents  of  all  input  cards  are  specified.  All 
possible  input  variables  are  listed  at  the  end  of  this  section  in  the 
order  of  appearance  in  the  input  deck.  The  input  format  for  all  cards  is 
shown  in  figure  E.3  and  the  item  numbers  also  refer  to  that  figure.  A 
sample  input  is  described  in  section  E,6  which  discusses  a sample  case. 

Item  1 

These  cards  provide  identification  of  the  run.  The  information  on 
them  is  printed  on  the  first  page  of  the  output. 

Item  2 

This  item  includes  the  variable  NFIN  which  specifies  the  configu- 
ration. The  following  configurations  may  be  studied: 

NFIN  = 1:  body-alone  configuration 

NFIN  = 2:  body-tail  configuration 

NFIN  = 3;  body-canard  configuration 

NFIN  = 4:  body-canard-tail  configuration 


r 


Item  3 


This  optional  item  is  read  if  NFIN  > 2 and  includes  integer 

control  parameters  for  the  first  set  of  fins.  The  first  six  specify 

whether  the  various  interference  parameters  are  to  be  input  by  the  user 

or  are  to  be  determined  by  the  program.  The  quantities  NCNC,  NCPXC, 

and  NCPyC  specify  the  number  of  entries  in  the  tables  which  describe  the 

wing-alone  characteristics.  Since  the  fins  are  assumed  to  be  unccimbered 

and  untwisted,  the  wing-alone  characteristics  should  be  input  for  positive 

a „ only . 
eq 

Also  in  this  item  is  NI,  the  number  of  integration  intervals  plus 
one,  used  in  the  Simpson's  rule  integration  routine  called  by  revflo. 

The  value  of  NI  must  be  odd.  A suitable  value  for  NI  depends  on  how 
close  a vortex  is  to  a fin.  If  a vortex  is  close  to  a fin,  NI  should  be 
at  least  51. 

The  quantity  NXOUTC  is  the  number  of  stations  for  which  vortex  posi- 
tions are  to  be  output.  A value  between  5 and  10  should  be  sufficient  to 
monitor  the  vortex  trajectories  over  the  canard  (first  finned)  section. 

The  quantity  NAFTC  is  used  if  NFKBC  and/or  NFXBC  are  zero  and 
> 1;  that  is,  and/or  Xg  are  to  be  determined  by  the  code.  For 
supersonic  speeds,  lift  on  the  body  due  to  the  presence  of  the  fins  can 
be  generated  downstream  of  the  fin  trailing  edges,  if  NAFTC  = 1,  the 
code  will  assume  that  there  is  sufficient  body  surface  downstream  of  the 
fin  trailing  edges  to  fully  develop  that  lift.  If  NAFTC  =2,  it  is 
assumed  that  there  is  no  body  after  the  fins,  if  NFIN  ^ 3,  always  use 
NAFTC  =1.  If  NFIN  = 2,  the  user  should  input  NAFTC  = 2 unless  it  is 
clear  that  significant  lift  carryover  will  be  developed  downstream  of 
the  fins. 

The  last  integer  variable  on  this  card  is  NVORT.  This  variable  is 
used  to  control  the  influence  of  the  nose  vortices.  It  has  been  observed 
(refs,  4 and  5)  that  nose  vorticity  for  < 20°  may  disperse  over  the 
canard  section.  Since  the  present  model  for  nose  vorticity  is  incap^le 
of  representing  such  a situation,  the  user  has  the  option  of  ignoring  the 
influence  of  nose  vorticity  downstream.  The  options  are: 

NVORT  “ 0:  nose  vortices,  if  present,  are  tracked  over  entire 

conf igurat ion 

NVORT  “ 1:  nose  vortices  are  ignored  downstream  of  the  leading 

edge  of  the  canard  root  chord 
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NVORT  = 2:  nose  vortices  are  ignored  downstream  of  the  trailing 

edge  of  the  canard  root  chord 

Item  4 

This  optional  item  is  read  if  NFIN  ^ 3.  The  quantity  NXOUTB  is 
equal  to  the  number  of  integration  intervals  used  in  the  calculations  of 
afterbody  loading.  It  is  also  the  number  of  axial  stations  at  which 
vortex  positions  and  strengths  are  output.  As  a rule  of  thumb,  use  a 
value  of  NXOUTS  greater  than  the  length  of  the  afterbody  in  body  diameters. 
If  the  vortex  positions  change  drastically  between  output  stations,  KXOUTB 
should  be  increased. 

Item  5 

This  optional  item  is  read  if  NFIN  = 4.  The  parameters  have  the 
same  function  as  those  in  Item  3.  For  excimple,  NXOUTT  has  the  same 
purpose  for  the  tail  section  as  NXOUTC  has  for  the  canard  section. 

Item  6 

The  logical  quantity  TURB  is  used  to  determine  which  branch  of  the 
crossflow  drag  coefficient  table  is  to  be  used  in  NOSE  and  AFTBOD.  For 
crossflow  Mach  numbers  below  0.6,  the  laminar  and  turbulent  values  of 
differ.  If  there  is  doubt  about  which  type  of  flow  separation  is 
present,  use  TURB  = .TRUE. 

The  logical  variables  NEARC  and  NEART  are  used  to  control  the 
influence  of  the  fins  on  vortex  trajectories  over  a finned  section.  If 
these  variables  are  set  equal  to  .TRUE.,  the  infloence  of  the  fins  is 
ignored.  This  option  has  given  the  best  agreement  with  experiment. 

Item  7 

This  item  contains  flight  conditions  and  fin  deflection  angles. 

Item  8 

This  optional  item  is  read  if  NFIN  >2.  It  is  possible  that  an 
indefinite  condition  may  result  if  PHI  is  zero  and  the  fins  are 
undeflected.  It  is  recommended  that  PHI  =0.01  be  used  for  an 
unbanked  configuration  with  undeflected  fins. 

Item  9 

This  item  contains  some  required  reference  information.  The 
variable  XMC  is  the  moment  center  measured  from  MS  0. 
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Item  10 


This  optional  information  for  the  first  set  of  fins  is  input  if 
NFIN  ^2.  The  variable  CPXRC  is  the  reference  length  which  was  used  to 
nondimens ionalize  the  chordwise  location  of  the  fin  center  of  pressure 
for  the  wing-alone  tables.  Typically,  this  length  is  equal  to  the  root 
chord  of  the  fin.  Similarly,  CPYRC  is  the  reference  length  for  the 
lateral  location  of  the  fin  center  of  pressure.  This  value  is  usually 
the  exposed  semispan  of  the  fin.  The  last  variable  is  SRCAN  which  is 
the  reference  area  used  for  the  fin  normal-force  table.  It  is  usually 
the  wing-alone  planform  area. 

Item  11 

This  optional  information  is  input  if  NFIN  = 4.  It  is  similar  to 
Item  10  and  is  the  corresponding  information  for  the  second  set  of  fins. 

Item  12 

This  information  is  input  if  NFIN  ^ 2.  It  is  concerned  with  control 
of  vortex  trajectory  calculations  and  with  the  influence  a vortex  may 
have  on  a fin.  The  first  variable,  RO,  enables  the  user  to  approximate 
the  effects  of  viscous  cores  for  the  vortices  in  determining  the  veloci- 
ties which  may  be  induced  on  the  fins,  if  RO  is  positive,  it  is  inter- 
preted as  the  ratio  of  the  vortex  core  radius  to  the  body  radius.  If 
RO  is  negative,  its  magnitude  is  interpreted  as  the  maximum  allowed 
fraction  of  the  free-strecim  velocity  which  can  be  induced  normal  to  the 
fins.  The  best  agreement  with  experiment  has  been  obtained  with 

RO  = -0.1  for  < 1 and  RO  = -0.4  for  M >1. 

00  00 

Subroutine  CRUTRJ  calculates  the  trajectories  of  the  vortices  in  the 
presence  of  fins.  It  uses  a differential  equation  solver,  DASCRU,  to 
solve  the  equations  of  motion  of  the  vortices.  The  input  variable,  DXI, 
is  the  initial  integration  interval  required  by  DASCRU  for  the  trajectory 
calculations.  It  is  nondimensionalized  by  the  body  radius.  The  initial 
interval  can  be  automatically  reduced  by  DASCRU  by  as  much  as  a factor 
of  100  in  an  attempt  to  provide  the  desired  accuracy.  A value  of  DXI 
between  0.02  and  0.05  should  work  for  most  cases.  When  two  vortices  get 
very  close  together,  or  when  a vortex  gets  very  near  the  wing  or  body, 
a smaller  value  of  DXI  may  be  required. 

The  desired  accuracy  in  the  vortex  trajectory  calculations,  referred 
to  in  the  previous  paragraph,  is  controlled  through  the  input  variable 
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EPS.  This  variable  also  controls  the  vortex  strength  calculations  of 
the  afterbody  vortices.  The  value  of  EPS  should  be  set  equal  to  one  half 
the  desired  accuracy.  For  instance,  if  EPS  is  set  to  5x10“^,  this  usually 
results  in  a solution  accurate  to  three  significant  figures. 

Experimental  observations  show  that  when  two  vortices  of  like  sign 
become  very  close  to  one  another  they  often  combine  or  coalesce  to  form 
one  vortex.  This  situation  is  included  in  the  program  through  the 
variable  RVORT.  If  vortices  are  not  combined  and  they  get  too  close, 
then  computation  time  can  increase  dramatically.  The  variable  RVORT  is 
entered  as  a separation  distance  in  terms  of  body  radii.  When  two  vor- 
tices of  like  sign  are  within  this  distance  of  each  other,  they  are 
combined  to  form  one  vortex  of  strength  equal  to  the  two  original 
strengths  and  at  a position  between  the  two  proportional  to  the  original 
strengths.  A value  of  0.5  for  RVORT  has  been  found  to  be  satisfactory. 

Item  13 

This  item  is  also  optional  and  contains  the  last  of  the  program 
control  information.  The  first  two  variables  in  this  item  govern  the 
shedding  of  the  afterbody  vortices.  Hence,  this  item  must  be  entered  if 
an  afterbody  is  present  behind  the  first  set  of  fins  (i.e.,  if  NFIN  ^ 3). 

The  first  variable  in  this  item  is  RVOA,  which  is  the  radial  distance 
from  the  body  centerline  in  body  radii  at  which  the  afterbody  vortices 
are  to  be  placed  initially.  The  suggested  range  of  values  for  RVOA  is 
1.2  to  1.5.  The  results  presented  in  this  report  were  computed  witn 
RVOA  = 1.2.  If  RVOA  is  too  small,  the  afterbody  vortices  tend  to  be 
captured  by  their  own  images. 

The  second  control  variable  in  this  item  is  OMEGA.  This  is  an 
angle,  measured  from  the  axis  within  which  the  crossflow  velocity 

vector  must  lie  in  order  for  afterbody  vortices  to  be  formed.  This 
crossflow  velocity  vector  is  determined  by  first  computing  the  upwash 
and  sidewash  at  the  body  eixis  due  to  the  external  vortices.  The  quantity 
sin  is  then  added  to  the  upwash  to  account  for  the  flow  of  the  free 
stream  around  the  cylindrical  body.  The  resultant  vector  is  the  crossflow 
velocity  vector.  If  OMEGA  is  set  equal  to  zero,  no  afterbody  vortices 
will  form  unless  the  flow  is  perfectly  symmetrical  about  the  plane  con- 
taining the  wind  and  body  axes.  If  OMEGA  is  set  equal  to  180°,  afterbody 
vortices  may  form  regardless  of  the  strength  and  direction  of  the  cross- 
flow  velocity.  A value  of  OMEGA  = 45°  is  recommended. 


The  last  variable  on  this  card  is  XTLE  which  is  the  axial  location 
of  the  end  of  the  afterbody  section. 

Item  14  1 I 

— 

This  item  describes  the  nose  section  back  to  the  first  set  of  fins. 

I 

The  first  card  contains  the  number  of  entries,  NNOSE,  in  the  table  of  nose  1 

coordinates  and  an  integer  variable,  NCA,  indicating  whether  the  slope  of  j 

the  linear  normal-force  curve,  is  to  be  read  in.  If  NCA  equals  zero,  ^ 

then  Cjj  is  assigned  the  value  2.0  by  the  program.  Otherwise,  a value 

ot  ) 

for  Cjj^  is  to  be  read  in  two  cards  later.  ^ 

Following  NNOSE  and  NCA,  the  two  variables  TIPRAD  and  ETAN  are 
entered.  The  variable  TIPRAD  is  the  radius  of  the  spherical  nose  tip. 

For  the  case  of  a pointed  nose,  this  value  is  zero.  If  the  nose  is  | 

pointed,  then  ETAN  is  the  nose  half-angle.  If  the  nose  is  blunted,  then  1 

ETAN  is  the  angle  between  the  body  axis  and  the  tangent  to  the  nose  at 
the  juncture  of  the  spherical  cap  and  the  rest  of  the  nose  (see  ref.  3, 

P.  11).  1 

If  NCA  is  not  equal  to  zero,  the  linear  normal-force  coefficient  ’■ 

slope  (DCNDA)  is  read  in  next.  Following  this  optional  input,  the  nose 
coordinates  are  entered,  first  the  axial  coordinates,  then  the  corres-  i 

ponding  radial  coordinates.  All  axial  locations  are  to  be  measured  from  | 

missile  station  zero  (MS  0).  Thus,  the  value  of  XNOSE(l)  is  the  cixial  ; 5 

location  of  the  nose  tip  measured  from  MS  0.  1 i 

j t 

Item  15  j I 

If  NFIN  = 1,  then  the  body  length  is  entered  in  this  item.  The  j ^ 

program  defines  the  nose  as  extending  to  the  leading  edge  of  the  root  I | 

chord  of  the  first  set  of  fins,  XCLE.  Thus,  for  the  body-alone  case,  | 1 

the  body  length  is  entered  as  XCLE  in  this  item.  g 

Item  16  I 

This  optional  item  is  entered  if  NFIN  > 2.  This  item  contains  I 

information  describing  the  first  set  of  fins.  On  the  first  card  of  this  a 

item  are  the  relevant  geometric  parameters  that  describe  the  fins.  If  | 

a fin  with  a rectangular  planform  is  to  be  specified,  the  user  must  set 
the  leading-edge  sweep  to  some  small  positive  value.  This  should  be 

accomplished  by  making  XCTIPL  > XCLE  by  a small  number,  say  0.01.  ^ 

Following  this  geometric  information  are  three  empirical  curves  that  j 

describe  the  fin  wing-alone  characteristics.  The  first  curve  is  the  t 
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i 


.J 


I 


! 

I 

1 

i normal-force  curve,  entered  as  corresponding  values  of  angle  of  attack 

I 

(ALFCNC)  and  normal-force  coefficient  (CNWC) . There  are  NCNC  pairs  of 
values.  The  second  empirical  curve  is  the  chordwise  center  of  pressure 
location  as  a function  of  angle  of  attack.  These  axial  locations  are 
measured  from  the  leading  edge  of  the  root  chord.  The  NCPXC  values  are 
entered  in  the  same  fashion  as  the  normal-force  curve.  Next,  the  table 
for  the  spanwise  location  of  the  wing-alone  center  of  pressure  is 
entered.  The  spanwise  distance  is  measured  from  the  root  chord.  There 
are  NCPYC  pairs  of  values. 

Following  these  empirical  curves,  any  or  all  of  the  interference 
factors  for  this  set  of  fins  may  be  entered.  As  discussed  previously, 
the  method  of  calculation  of  each  interference  factor  is  controlled  by 
an  integer  flag.  if  the  value  of  the  flag  is  zero,  then  the  corresponding 
interference  factor  is  not  entered  but  is  calculated  internally.  If  the 
value  of  the  flag  is  one,  then  the  interference  parameter  is  entered  here 
for  the  first  set  of  fins  and  at  the  end  of  Item  17  for  the  second  set 
of  fins.  The  integer  flags  corresponding  to  tlie  interference  factors 
for  the  set  of  fins  of  Item  16  were  entered  in  Item  3.  Each  card  is  read 
if  any  of  the  information  on  it  is  needed.  The  quantities  which  c*re  not 
needed  are  ignored. 

Item  17 

This  item  is  similar  to  Item  16  in  that  it  contains  the  corresponding 

1 

information  for  the  second  set  of  fins.  Hence,  this  item  is  entered  only 
when  NFIN  = 4.  In  addition  to  geometrical  information,  the  first  card 
contains  the  interdigitation  angle.  This  is  the  angle  by  which  the  tail 
fins  are  rolled  with  respect  to  the  canard  fins.  The  axial  location  of 
the  tail  fin  root  chord  leading  edge  was  entered  in  Item  13.  Just  as  in 
Item  16,  the  three  empirical  curves  describing  the  wing-alone  characte- 
ristics are  entered  next.  The  last  two  cards,  if  they  are  entered, 
contain  the  fin  interference  factors.  The  integer  flags  that  determine 
whether  or  not  the  interference  factors  for  this  fin  are  entered  are  in 
Item  5.  Each  card  is  read  if  any  of  the  information  on  it  is  needed. 

The  quantities  which  are  not  needed  are  ignored. 

Item  18 

This  card  ends  the  process  of  entering  data.  It  should  be  the  last 
card  and  follow  all  the  data  cards  for  the  case(s)  to  be  run.  The 
computer  program  stops  the  search  for  more  data  and  the  run  is  finished. 
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Item  No.  1: 


(1)  Format  (1615),  1 card 
Column  Number  S \ 

Pro9ram  Variable  NCARDS  \ 


(2)  Format  (20A4) , NCARDS  cards 
Column  Number  |l-80 


Program  Variable 


Item  No.  2: 


Column  Number 


Item  No.  3:  Optional 


5 

10 

NFIN 

NALFA 

Column  Number 


(1)  Format  (1615) 

I sf  l3 


1C 

15 

20 

25 

30 

NFLKW 

»IFKABC 

NFKD 

NFKBC 

NFXBC 

50 

55 

60 

65 

70 

NI 

NPHI 

dXOOTC 

NAFTC 

NVORT 

Item  No.  4:  Optional 


Column  Number 

Program  Variable  NXOOTB 


(1)  Format  (1615) 

I 


Figure  B.3(a).  Input  formats  for  computer  program  MLOAOS. 


Item  No.  5:  Optional 


Column  Number 
Program  Variable 


(1)  Format  (1615) 


5 

10 

15 

20 

25 

30 

35 

NFKWT 

NFKABT 

NFKBT 

NFXBT 

NCNT 

NCPXT 

NCPYT 

40 

45 

I 

NXOUTT 

NAFTT 

] 

Item  No.  6: 

(1)  Format  (16L5) 


Column  Number 

5 

10 

15 

3 

Program  Varietble 

TURB 

NEARC 

NEART 

□ 

Item  No.  7: 

(1)  Format  (8F10.5) , 1 card 


Column  Number 

10 

20 

30 

40 

50 

7 

Program  Variedjle 

FMACH 

DELTA  (1) 

DELTA (2) 

DELTA (3) 

DELTA (4) 

7 

(2)  Format  (8F10.5) , 8 values  of  ALFAC  per  card 


Column  Number 

10 

20 

10  X NALFA 

1 

Program  VarieU^le 

ALFAC (1) 

ALFAC  (2) 

. . . 

ALFAC  (NALFA 

J 

Item  No.  8:  Optional 

(1)  Format  (8F10.5) , 8 values  of  PHI  per  card 


Column  Number 

10 

20 

• • • 

10  X NPHI 

7 

Program  VarieOBle 

PHI(l) 

PHI  (2) 

. . . 

PHI (NPHI ) 

Figure  E.3(b).  Continued. 
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Colunm  Humber 
Program  Variable 


Column  Number 
Program  variable 


Column  Number 
Program  Variable 

Item  No.  18: 

Column  Number 
Program  Variable 


(6)  Format 

(8F10.5)  , 

8 values  of  ACPYT  per 

m 

10 

. . . 

lOxNCPVT 

1 

ACPYT(l) 

ACPVT(2) 

. . . 

ACPYT (NCPYT) 

d 

10 

. . . 

10  X NCPYT 

CPYT(l) 

CPYT(2) 

. . . 

CPYT  (NCPYT) 

(8)  Optional,  Format  (8F10.5),  1 card 


10 

20 

30 

40 

50 

60 

70 

KWT 

KBT 

XBT 

KABT(l) 

KABT(2) 

KABT(3) 

KABT(4) 

(1)  Format 
5 

999 

(1615) 

K 

Figure  E.3(g).  Concluded 
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E.4.  SYMBOL  LISTING 

This  section  lists  all  input  varialles  in  the  order  in  which  they 
are  entered,  gives  the  algebraic  symbol  if  one  exists,  and  a brief 
definition.  Refer  to  figure  E.3  for  the  proper  input  format. 


Item  2 
NFIN 


NALFA 

Item  3 (optional) 
NFKWC 

NFLKW 

NFKABC 

NFKD 


Comments 

Alphanumeric  information  to  identify  the 
run. 

Number  of  cards  used  to  identify  the  run; 
NCARDS  ^ 1. 

(NCARDS)  cards  of  alphanumeric  information 
for  identification  of  the  run; 

1 < I < NCARDS. 


Integer  specifying  type  of  configuration. 
NFIN  = 1 Body- alone. 

NFIN  = 2 Body-tail. 

NFIN  = 3 Body-canard. 

NFIN  = 4 Body-canard-tail. 

Number  of  angles  of  attack  for  which  cal- 
culations are  to  be  made;  1 < NALFA  <10. 

Optional  input  to  be  entered  if  NFIN  ^ 2. 

Integer  flag  indicating  how  for 

canard  section  is  to  be  determined. 

NFKWC  = 0 Determined  by  code. 

NFKWC  =1  Kw  is  input. 

Integer  flag  indicating  how  for 

canard  section  is  to  be  determined. 

NFLKW  = 0 Determined  by  code. 

NFLKW  =1  ^ input. 

Integer  flag  indicating  how  for 

canard  section  is  to  be  determined. 

NFKABC  = 0 Determined  by  code. 

NFKABC  = 1 K^p  is  input. 

Integer  flag  indicating  how  for 

canard  section  is  to  be  determined. 

NFKD  = 0 Determined  by  code. 

NFKD  =1  K ^ is  input. 


Program 

variable 


Algebraic  Symbol 
(If  Applicable) 


Comments 


NFKBC 

NFXBC 

NCNC 

NCPXC 

NCPyC 

NI 

NPHI 

NXOUTC 

NAFTC 

ISVORT 


Item  4 (optional) 


Integer  flag  indicating  how  Kg  for 
canard  section  is  to  be  determined. 

NFKBC  = 0 Determined  by  code. 

NFKBC  = 1 Kg  is  input. 

Integer  flag  indicating  how  Xg  for 
canard  section  is  to  be  determined. 

NFXBC  = 0 Determined  by  code. 

NFXBC  =1  Xg  is  input. 

Number  of  entries  in  table  for  canard 
wing-alone  normal  force  coefficient; 

2 ^ NCNC  ^ 20. 

Nuniber  of  entries  in  table  for  the  chord- 
wise  location  of  the  canard  wing-alone 
center  of  pressure;  2 < NCPXC  <20. 

Number  of  entries  in  table  for  the  span- 
wise  location  of  the  canard  wing-alone 
center  of  pressure;  2 < NCPYC  20. 

One  plus  the  number  of  intervals  to  be 
used  in  the  Simpson's  rule  integration 
package  in  REVFLO;  must  be  odd;  1^NI<;^  99. 

Number  of  roll  angles  for  which  calcula- 
tions are  to  be  made;  1 < NPHI  <10. 

Number  of  axial  stations  along  the  canard 
at  which  vortex  locations  are  to  be 
output . 

Integer  flag  indicating  whether  afterbody 
is  present  downstream  of  the  first  set 
of  fins  for  computation  of  body  carryover 
lift;  ignored  if  NFKBC  and  NFXBC  ^ 0. 
NAFTC  = 1 Afterbody  present. 

NAFTC  = 2 NO  afterbody  present. 

Integer  flag  indicating  how  far  along  body 
influence  of  nose  vortices  is  to  be  felt. 
NVORT  = 0 Influence  of  nose  vortices 
felt  along  entire  body. 

NVORT  = 1 Influence  of  nose  vortices 
felt  up  to  leading  edge  of 
canard  root  chord. 

NVORT  = 2 influence  of  nose  vortices 

felt  up  to  trailing  edge  of 
canard  root  chord. 

Optional  input  to  be  read  in  if  NFIN  > 3. 


NXOUTB 


Number  of  axial  stations  along  the  after- 
body at  which  vortex  locations  are  to  be 
output  (see  Item  4 of  section  E.3) . 


Program 

variable 


Algebraic  Symbol 
(If  Applicable) 


Comments 


Item  5 (optional) 
NFKWT 

NFKABT 

NFK3T 

NFXBT 

NCNT 

NCPXT 

NCPYT 

NXOUTT 

NAFTT 


Item  6 
TURB 


Optional  input  to  be  read  in  if  NFIN  = 4. 

Integer  flag  indicating  how  for  tail 

section  is  to  be  determined. 

NFKWT  = 0 Determined  by  code. 

NFKWT  =1  Kw  input. 

Integer  flag  indicating  how  K^p  for  tail 
section  is  to  be  determined. 

NFKABT  = 0 Determined  by  code. 

NFKABT  = 1 K^p  is  input. 

Integer  flag  indicating  how  Kg  for  tail 
section  is  to  be  determined. 

NFKBT  = 0 Determined  by  code. 

NFKBT  =1  Kg  is  input. 

Integer  flag  indicating  how  Xg  for  tail 
section  is  to  be  determined. 

NFXBT  = 0 Determined  by  code. 

NFXBT  =1  Xg  is  input. 

Number  of  entries  in  table  for  tail  wing- 
alone  normal  force;  2 ^ NCNT  ^ 20. 

Number  of  entries  in  table  for  the  chord- 
wise  location  of  the  tail  wing-alone 
center  of  pressure;  2 < NCPXT  ^ 20. 

Number  of  entries  in  table  for  the  span- 
wise  location  of  the  tail  wing-alone 
center  of  pressure;  2 < NCPYT  < 20. 

Number  of  axial  stations  along  the  tail 
at  which  vortex  locations  are  to  be 
output. 

Integer  flag  indicating  whether  afterbody 
is  present  downstream  of  second  set  of 
fins  for  calculation  of  body  carryover 
lift;  ignored  if  NFKBT  and  NFXBT  ^ 0. 
NAFTT  = 1 Afterbody  present. 

NAFTT  =2  No  afterbody  present. 


Logical  variable  stating  whether  crossflow 
on  body  is  laminar  or  turbulent. 

TURB  = .TRUE.  Crossflow  is  turbulent. 

TURB  = .FALSE.  Crossflow  is  laminar. 
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Program  Algebraic  Symbol 

Variable  (If  Applicable)  Comments 


NEARC 

Logical  variable  stating  whether  canard 
fins  are  to  be  ignored  for  vortex 
tracking  over  canard  section;  NEARC 
is  ignored  if  NFIN  < 2. 

NEARC  = .TRUE.  Ignore  fins. 

NEARC  = .FALSE.  Include  influence  of 

fins. 

NEART 

Logical  variable  stating  whether  tail  fins 
are  to  be  ignored  for  vortex  tracking 
over  tail  section;  NEART  is  ignored  if 
NFIN  < 4. 

NEART  = .TRUE.  Ignore  fins. 

NEART  = .FALSE.  Include  influence  of 

fins. 

Item  7 

FMACH 

Free-stream  Mach  nvmiber. 

DELTA (J) 

6 (J) 

Deflection  angles  in  degrees  of  canard  fins; 
1 < J < 4;  J = number  of  canard  fin. 

ALFAC (K) 

ac(K) 

Body  angle  of  attack  in  degrees; 
1 < K < NALFA. 

Item  8 (optional) 

Optional  input  to  be  read  in  if  NFIN  ^2. 

PHI  (L) 

0(L) 

Bank  angle  in  degrees;  angle  between  z 
and  2q  axes;  positive  measured  clockwise 
viewed  from  rear;  1 < L < NPHI. 

Item  9 

LROUT 

^r 

Reference  length  used  in  center  of  pressure 
calculations . 

SROUT 

Sr 

Reference  area  used  in  force  calculations. 

XMC 

Moment  center  of  missile,  dimensional. 

A 

a 

Radius  of  missile,  dimensional. 

Item  10  (optional) 

Optional  input  to  be  read  in  if  NFIN  ^ 2. 

CPXRC 

Reference  length  used  in  table  for  chord- 
wise  location  of  canard  wing-alone  center 
of  pressure,  dimensional. 

CPYRC 

Reference  length  used  in  table  for  span- 
wise  location  of  canard  wing-alone  center 
of  pressure,  dimensional. 

SRCAN 

Reference  area  used  in  table  for  canard 

wing-alone  normal-force  coefficient, 
dimensional . 
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Program 

variable 


Algebraic  Symbol 
(If  Applicable) 


Comments 


Item  11  (optional) 
CPXRT 

CPYRT 

SRTAIL 

Item  12  (optional) 
RO 


DXI 

EPS 

RVORT 

Item  13  (optional) 

RVOA  r^a 

OMEGA 

XTLE 

Item  14 
NNOSE 


Optional  input  to  be  read  in  if  NFIN  = 4. 

Reference  length  used  in  table  for  chord- 
wise  location  of  tail  wing-alone  center 
of  pressure,  dimensional. 

Reference  length  used  in  table  for  span- 
wise  location  of  tail  wing-alone  center 
of  pressure,  dimensional. 

Reference  area  used  in  table  for  tail 
wing-alone  normal-force  coefficient, 
dimensional . 

Optional  input  to  be  read  in  if  NFIN  ^ 2. 

If  rO  ^ 0,  RO  is  the  ratio  of  the  vortex 
core  radius  to  the  body  radius;  if 
RO  ^ 0,  RO  is  the  negative  of  the  maxi- 
mum fraction  of  the  free-stream  velocity 
which  may  be  induced  normal  to  the  fins. 

Mciximum  value  of  integration  interval  in 
body  radii  for  vortex  tracking. 

Accuracy  criterion  used  in  DASCRU. 

If  the  separation  distance  in  body  radii 
between  vortices  of  like  sign  is  less 
than  RVORT,  vortices  are  combined 
(afterbody  section  only) . 

Optional  input  to  be  read  in  if  NFIN  ^ 3. 

Distance  from  body  centerline  at  which 
afterbody  vortices  are  positioned 
initially,  nondimens ionalized  by  body 
radius. 

Angle  measured  from  Zq  axis  within  which 
the  crossflow  velocity  vector  must  lie  in 
order  for  afterbody  vortices  to  be  formed, 
degrees . 

Location  of  end  of  afterbody  section; 
position  of  leading  edge  of  tail  root 
chord  if  second  set  of  fins  are  present 
and  axial  position  of  body  base  if  only 
one  set  of  fins  present. 


Number  of  entries  in  the  table  of  nose 
coordinates;  2 < NNOSE  ^ 20. 
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Program 

Variable 


Algebraic  Symbol 
(If  Applicable) 


Comments 


NCA 


TIPRAD 

ETAN  T) 


DCNDA  Cnq 

XNOSE  (M) 


RNOSE(M) 


Item  15  (optional) 
XCLE 


Item  16  (optional) 


ARC  iR 

SPANC 

m 


XCHL 

XCLE 


XCTIPL 


XCTE 

XFC 


Integer  flag  specifying  whether  of 

nose  is  to  be  entered. 

NCA  = 0 not  entered. 

NCA  / 0 entered. 

Nose  tip  radius,  dimensional. 

Half  angle  of  body  nose  for  pointed  body; 
or  angle  between  tangent  to  nose  at 
juncture  of  spherical  cap  and  rest  of 
nose  and  the  body  axis,  degrees. 

Slope  of  nose  normal-force  coefficient  at 
= 0;  entered  if  NCA  / 0. 

Axial  location  entries  in  nose  coordinate 
table  from  a nose  tip;  1 ^ M < NNOSE; 
XNOSE (1)  is  axial  location  of  nose  tip 
from  MS  0,  dimensional. 

Corresponding  radial  location  entries  in 
nose  coordinate  table;  1 < M ^ NNOSE, 
dimensional . 

Optional  input  to  be  read  in  if  NFIN  = 1. 

Length  of  body,  dimensional. 

Optional  input  to  be  read  in  if  NFIN  ^ 2. 

Aspect  ratio  of  canard  fins. 

Maximum  semi span,  measured  from  body 
centerline,  of  canard  fins,  dimensional. 

Axial  distance  to  canard  hinge  line, 
measured  from  MS  0,  dimensional. 

Axial  distance  to  leading  edge  of  canard 
root  chord,  measured  from  MS  0, 
dimensional . 

Axial  distance  to  leading  edge  of  canard 
tip  chord,  measured  from  MS  0, 
dimensional;  XCTIPL  > XCLE  (see  Item  16 
of  section  0.3) . 

Axial  distance  to  canard  trailing  edge, 
measured  from  MS  0,  dimensional. 

Axial  distance  to  area  centroid  of  canard 
fin  planform,  measured  from  MS  0, 
dimensional . 
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Program 

Variable 


Algebraic  Symbol 
(If  Applicable) 


Comments 


ALFCNC (N) 
CNWC (N) 
ACPXC(II) 
CPXC(II) 


ACPyC  (IJ) 
CPYC  ( IJ) 


KWC 


Angle  of  attack  entries  in  canard  wing- 
alone  normal-force  table,  degrees; 

2 < N < NCNC . 

Corresponding  normal-force  coefficient 
entries  in  canard  wing-alone  normal-force 
table;  2 < N ^ NCNC. 

Angle  of  attack  entries  in  canard  wing- 
alone  table  for  chordwise  location  of 
center  of  pressure,  degrees,  2 ^ II  ^ NCPXC. 

Corresponding  entries  for  chordwise  locations 
of  center  of  pressure  in  canard  wing-alone 
table,  nondimensional ; 2^11^  NCPXC, 
measured  from  leading  edge  of  canard  root 
chord . 

Angle  of  attack  entries  in  canard  wing- 
alone  table  for  spanwise  location  of  center 
of  pressure,  degrees;  2 < IJ  ^ NCPYC. 

Corresponding  entries  for  spanwise  locations 
of  center  of  pressure  in  canard  wing-alone 
table,  nondimensional;  2 ^ IJ  ^ NCPYC, 
measured  from  fin  root  chord. 

Wing-body  interference  factor  for  canard 
normal  force . 


LKW 

Wing-body  interference  factor  for  canard 
fin  deflection. 

KBC 

Wing-body  interference  factor  for  lift 
carryover  from  canard  fins  to  body. 

XBC 

Center  of  pressure  for  Kg  load  based  on 
fraction  of  root  chord  measured  from 
leading  edge  of  canard  root  chord. 

KABC  (IK) 

V7ing-body  interference  factor  associated 
with  a3  coupling  on  canard  fins; 

1 ^ IK  ^ 4;  IK  = number  of  canard  fin. 

KD(IL) 

K^(IL) 

Wing -body  interference  factor  associated 
with  a6y  and  coupling  on  canard 

fins;  1 < IL  ^ 4;  IL  = number  of  canard 
fin. 

Item  17  (optional)  Optional  input  to  be  read  in  if  NFIN  = 4. 


ART 

Aspect  ratio  of  tail  fins. 

SPANT 

®m 

Maximum  semispan,  measured 
centerline,  of  tail  fins. 

from  body 
dimensional . 

178 


Program  Algebraic  Symbol 
variable  (If  Applicable) 


Comments 


XTTIPL 

XTTE 

XFT 

PHIT  0^ 

ALFCNT(IM) 

CNWT ( IM) 

ACPXT(IN) 

CPXT(IN) 

ACPYT(KI) 

CPyT(KI) 


KWT 

KBT 

XBT 

KABT(KK) 

Item  18 

999 


Axial  distance  to  leading  edge  of  tail  tip 
chord,  measured  from  MS  0,  dimensional; 
XTTIPL  > XTLE  (see  Item  16  of  section  E.3) . 

Axial  distance  to  tail  fin  trailing  edge, 
measured  from  MS  0,  dimensional. 

Axial  distance  to  area  centroid  of  tail  fin 
planform,  measured  from  MS  0,  dimensional. 

Interdigitation  angle  between  tail  fins  and 
canard  fins,  measured  clockwise  from  canard 
fin  no.  1 to  tail  fin  no.  1,  degrees 

Angle  of  attack  entries  in  tail  wing-alone 
normal-force  table,  degrees;  2 ^ IM  < NCNT. 

Corresponding  normal-force  coefficient 
entries  in  tail  wing-alone  normal-force 
table;  2 ^ IM  ^ NCNT. 

Angle  of  attack  entries  in  tail  wing-alone 
table  for  chordwise  location  of  center  of 
pressure,  degrees;  2 < IN  ^ NCPXT,  measured 
from  leading  edge  of  tail  root  chord. 

Corresponding  entries  for  chordwise  loca- 
tions of  center  of  pressure  in  tail  wing- 
alone  table,  nondimensional;  2 ^ IN  ^ NCPXT. 

Angle  of  attack  entries  in  tail  wing-alone 
table  for  spanwise  location  of  center  of 
pressure,  degrees;  2 ^ KI  ^ NCPYT. 

Corresponding  entries  for  spanwise  loca- 
tions of  center  of  pressure  in  tail  wing- 
alone  table,  nondimensional;  2 ^ KI  ^ 

NCPYT,  measured  from  fin  root  chord. 

Wing-body  interference  factor  for  tail 
normal  force. 

Wing-body  interference  factor  for  lift 
carryover  from  tail  fin  to  body. 

Center  of  pressure  for  Kg  based  on  per- 
centage of  root  chord  measured  from 
leading  edge  of  tail  root  chord. 

Wing-body  interference  factor  associated 
with  ap  coupling  on  tail  fins; 

1 ^ KK  ^ 4;  KK  * number  of  tail  fin. 


This  card  causes  the  program  to  stop 
searching  for  more  data  and  the  run  is 
stopped . 
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E.5.  DESCRIPTION  OF  OUTPUT 

This  section  describes  the  output  of  the  computer  program.  The 
contcmts  of  each  of  the  output  items  is  specified  and  discussed,  A 
sample  output  is  included  in  the  next  section  which  discusses  the 
sample  case. 

The  output  consists  of  a minimum  of  four  items  and  a maximum  of 
seven,  depending  upon  the  type  of  configuration.  The  first  item  consists 
of  the  input  data  and  takes  up  the  first  three  pages.  The  run  identifi- 
cation and  missile  geometry  appear  on  the  first  page.  The  second  page 
contains  reference  lengths,  reference  areas,  and  flight  conditions. 
Included  with  the  flight  conditions  arc  the  values  of  TURB,  NEARC,  and 
NEART  listed  under  the  headings,  "turbulent  flow,"  "vortices  are  near 
canard,"  and  "vortices  are  near  tail,"  respectively.  The  third  and  last 
page  of  this  item  contains  the  tables  for  the  wing-alone  characteristics 
for  both  sets  of  fins.  The  program  operation  variables  are  also  given 
on  this  page. 

The  second  output  item  appears  on  page  four  and  is  a listing  of  all 
interference  parameters  used  in  the  program  and  their  source.  The  number 
xero  for  the  source  of  an  interference  parameter  indicates  that  that 
particular  parameter  was  calculated  by  the  code.  If  the  value  for  the 
source  is  one,  the  interference  parameter  was  input  by  the  user. 

Page  five  marks  the  beginning  of  the  force  and  moment  results. 

Three  items  of  information  are  given  at  the  beginning  of  the  output  for 
each  of  the  individual  sections;  (1)  the  name  of  the  particular  section 
(e.g.,  nose);  (2)  the  value  of  and  (3)  the  bank  angle,  0.  Following 

the  output  for  the  last  section  of  the  configuration  is  a summary  giving 
total  forces  and  moments.  The  forces  and  moments  are  all  presented  in 
coefficient  form  with  SROUT  as  the  reference  area  and  LROUT  as  the  refer- 
ence length.  In  addition,  the  force  and  moment  coefficients  are  listed 
as  components  along  the  appropriate  Cartesian  axes.  The  coefficients 
in  unrolled  coordinates  are  listed  as  C20,  normal  force;  CYO,  side  force, 
exo,  axial  force;  CMZO,  yawing  moment;  CMYO,  pitching  moment;  and  CMXO, 
rolling  moment.  The  same  convention  is  used  for  the  values  in  rolled 
coordinates.  Note  that  CZO  = CYO  = Cy,  CXO  = C^,  CMXO  = -C^, 

CMYO  = C|j^,  and  CMZO  = -C^-  Also  printed  is  the  contribution  to  the  lift 
and  drag  coefficient  from  each  section. 
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The  loads  on  the  nose  section  are  reported  on  the  fifth  page  of  the 
output  beginning  with  the  axial  starting  position  of  the  nose  vortices, 
if  any.  If  nose  vortices  have  formed,  their  nondimens ional  strength 
(normalized  by  2TrV^a)  and  their  positions  in  the  crossflow  plane  at  the 
end  of  the  nose  section  are  presented  next.  Then  the  nose  forces  and 
moments  are  printed. 

The  results  from  the  canard  section  (or  first  finned  section)  start 
on  page  six.  Presented  first  are  results  for  the  individual  fins  in  the 
absence  of  vortices.  For  each  fin  the  following  information  is  given: 

(1)  the  equivalent  angle  of  attack;  (2)  the  fin  normal-force  coefficient; 
and  (3)  the  location  of  the  fin  center  of  pressure  normalized  by  LROUT. 

The  value  CPX  is  measured  from  the  root  chord  heading  edge,  and  CEY  is 
measured  from  the  body  centerline. 

If  nose  vortices  have  formed  and  if  their  influence  over  the  canard 
section  is  to  be  calculated  (i.e.,  if  NVORT  is  not  equal  to  one),  the 
results  from  CRUTRJ  giving  the  trajectory  of  the  vortices  over  the  fin 
section  are  presented  next.  Positions  in  the  crossflow  plane  in  both 
rolled  and  unrolled  coordinates  are  given  at  NXOUTC+2  axial  locations, 
the  first  one  being  the  root  chord  leading  edge.  Other  output  locations 
are  at  intervals  of  (XCTE  - XCLE)/NXOUTC  with  the  positions  at  the 
axial  location  of  the  area  centroid  also  given.  It  is  the  vortex  posi- 
tions at  the  fin  area  centroid  that  are  used  in  computing  the  vortex 
induced  effects.  Immediately  following  the  vortex  trajectory  information 
are  the  vortex  induced  equivalent  angles  of  attack  and  corresponding 
spanwise  locations  of  the  centers  of  pressure  for  each  fin.  The  spanwise 
locations  of  the  centers  of  pressure  are  measured  from  the  body  centerline 
and  are  nondimens ional ized  by  the  body  radius. 

Next  are  listed  the  strengths  and  positions  of  the  canard  fin 
trailing  vortices  in  the  crossflow  plane  at  the  canard  fin  trailing  edge. 
One  vortex  is  generated  per  fin.  Total  individual  fin  loads  are  pre- 
sented next.  For  each  fin  the  total  equivalent  angle  of  attack  is  used 
to  determine  the  force  normal  to  the  fin  (CN) , the  location  of  the 
normal-force  center  of  pressure  (CPX  and  CPY) , the  fin  rolling  moment 
(CRM) , bending  moment  about  the  root  chord  (CBM) , and  hinge  moment  (CHM) . 
These  individual  fin  loads  are  then  sxunmed  and  added  to  the  forces  and 
moments  acting  on  the  body  in  the  presence  of  the  fins.  The  totals  are 
given  in  both  rolled  and  unrolled  coordinates.  If  a body-tail  missile 


is  being  run,  the  results  from  the  tail  section,  being  the  first  set  of 
fins,  appear  here. 

Following  the  canard  section  results  are  the  load  contributions  from 
the  afterbody  section.  These  results  are  presented  in  the  form  of  vortex, 
flow  field,  and  force  information  for  each  of  the  NXOUTB  axial  stations 
specified.  From  left  to  right  the  printout  reads:  X/A,  the  axial 

station,  nondimensionalized  by  the  body  radius;  I,  the  identification 
number  of  the  vortices;  GAMMA(I),  the  nondimens ional  strength  of  the 
vortices  at  X/A;  Y(I)/A,  the  nondimens ional  y^-coordinate  of  the  vortices; 
Z(I)/A,  the  nondimens ional  z^-coordinate  of  the  vortices;  FLOANG,  the 
inverse  sine  in  degrees  of  the  magnitude  of  the  crossflow  velocity  vector 
nondimensionalized  by  V^;  THETA,  the  polar  angle  in  degrees  of  the 
crossflow  velocity  vector  measured  from  the  positive  y^-axis  counter- 
clockwise; NORMAL  FORCE,  the  increment  of  normal  force  acting  on  the 
afterbody  over  the  preceding  interval;  and  SIDE  FORCE,  the  increment  of 
side  force  acting  on  the  afterbody  over  the  preceding  interval.  If  two 
vortices  are  combined  (i.e.,  if  two  vortices  of  like  sign  come  within  a 
distance  of  RVORT  of  each  other) , then  their  strengths  are  added  and  the 
indexing  of  the  vortices  is  rearranged  as  described  in  the  output.  The 
contribution  of  the  afterbody  section  to  the  total  loads  in  the  unrolled 
body  coordinate  system  is  presented  after  the  trajectory  information. 

Results  for  the  tail  section  appear  next.  The  format  is  the  same 
as  for  the  canard  section  with  the  exception  of  hinge  moment.  The 
contribution  of  the  tail  section  to  the  axial  force  is  always  zero 
because  the  tail  panels  do  not  deflect. 

The  final  item  of  output  is  the  summary  of  the  total  forces  and 
moments  on  the  configuration  in  both  rolled  and  unrolled  coordinates. 
Included  in  this  section  is  the  axial  center  of  pressure  due  to  the 
total  normal  force  (CEX)  and  the  axial  center  of  pressure  due  to  the 
total  side  force  (CPY) . 


E.6.  SAMPLE  CASE 


In  this  section  we  describe  a sample  case  to  illustrate  the  use  of 
the  computer  program.  The  case  consists  of  a body-canard-tail  configu- 
ration at  supersonic  speed,  unrolled,  at  an  angle  of  attack,  and  with 
yaw  control. 
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Figure  E.4.  Input  for  sample  case; 
body-canard-tail  configuration. 
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Figure  E. 5(c).  Continued 
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Figure  E.5(d).  Continued, 
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Figure  E.  5(e)  Continued, 
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The  configuration  is  shown  in  figure  3 (h)  of  the  main  text  and  is 
the  MICOM-NWC  model  of  reference  4 with  the  N„  nose,  the  C canard 
fins,  and  the  T tail  fins.  The  angle  of  attack  is  16.6°,  the  model 
is  at  a bank  angle  of  0°,  and  the  yaw  panels  (canard  panels  1 and  3)  are 
deflected  15°.  The  pitch  panels  are  not  deflected.  Mach  number  is  1.75 
and  the  crossflow  is  considered  turbulent.  The  values  of  five  interfer- 
ence factors  (KBC,XBC,KWT,KBT,XBT)  were  input.  These  values  were 
obtained  from  running  the  wing-body  computer  program  CRFWBD  described 
in  Appendix  F and  reference  3.  The  influence  of  nose  vortices  is  calcu- 
lated to  the  canard  trailing  edge  (NVORT  = 2) . Both  NEARC  and  NEART 
are  set  equal  to  .TRUE.  The  maximum  vortex  induced  velocity  on  the  fins 
is  limited  to  0.4  (RO  = -0.4) . Afterbody  vortices  (if  they  form) 
are  initially  positioned  at  1.2  radii  from  the  body  cixis  (RVOA  = 1.2). 
Finally,  vortices  of  like  sign  are  combined  along  the  afterbody  if  they 
come  within  a distance  of  0.5  body  radii  of  one  another  (RVORT  = 0.5) . 
The  input  deck  for  this  case  is  shown  in  figure  E.4.  The  output  is 
presented  in  figure  E.5. 


E.7.  PROGRAM  LISTING 

The  program  is  written  in  FORTRAN  IV  computer  language  (029  punch) . 
The  program  consists  of  a main  program,  MLOADS,  and  twenty-eight  sub- 
programs. Each  source  deck  is  identified  in  columns  73  through  80  by 
a three-character  identification  and  a three-digit  number  sequencing 
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APPENDIX  F ‘ 


MODIFICATIONS  TO  CRFWBD,  A COMPUTER 
PROGRAM  FOR  CALCULATING  AERODYNAMIC 
CHARACTERISTICS  OF  CRUCIFORM  WING-BODY 
COMBINATIONS  IN  SUPERSONIC  FLOW 


APPENDIX  F 

MODIFICATIONS  TO  CRFWBD,  A COMPUTER 
PROGRAM  FOR  CALCULATING  AERODYNAMIC 
CHARACTERISTICS  OF  CRUCIFORM  WING-BODY 
COMBINATIONS  IN  SUPERSONIC  FLOW* 


F.l.  INTRODUCTION 


This  appendix  describes  several  changes  and  extensions  to  the 
computer  progrcun  CRFWBD  first  discussed  in  Appendix  A of  reference  3. 

In  addition,  a study  was  made  to  determine  the  effect  on  fin  loading 
of  various  panel  sizes  used  to  cover  the  fin  planform  and  body  interfer- 
ence shell.  The  results  of  this  investigation  are  also  included.  Only 
the  modifications  to  the  program  CRFWBD  which  were  made  during  the  course 
of  the  present  study  will  be  discussed  in  this  appendix.  The  user  should 
consult  Appendix  A of  reference  3 for  a complete  description  of  the 
original  code.  The  modified  program  is  written  in  the  FORTRAN  IV  language 
for  the  CDC  7600  machine.  The  modified  routines  are  punched  in  029.  The 
user  should  note  that  the  original  routines  are  all  punched  in  026. 


F.2.  MODIFICATIONS  TO  CRFWBD 


The  first  change  to  the  program  was  to  remove  all  calls  to  the 
previously  system-supplied  subroutine  REQFL.  This  subroutine  computed 
the  actual  dimension  requirement  for  the  aerodynamic  coefficierit  matrix 
FVN.  In  the  current  version  this  matrix  is  given  the  fixed  dimension 
corresponding  to  the  maximum  number  of  constant  u-velocity  panels 
allowed  and  is  placed  in  large  core  memory.  The  maximum  number  of 
constant  u-velocity  panels  available  has  been  increased  from  250  to  350. 
Thus,  the  dimension  of  FVN  is  350x350  or  122,500. 

Several  modifications  were  made  to  subroutine  LOADS  to  compute  the 
lift  carryover  onto  the  body.  In  addition,  the  variable  LIL  (for  body 
i.nterference  length)  has  been  added  to  namelist  INPUT  to  allow  the  user 
to  select  the  length  (measured  from  the  leading  edge  of  the  fin  root 
chord)  over  which  the  wing  exerts  lift  carryover  onto  the  body.  In 
general,  the  minimum  length  that  should  be  chosen  is  determined  by  the 
intersection  of  the  Mach  cone  from  the  trailing  edge  of  the  root  chord 
with  the  body  centerline, 

*The  work  described  in  this  appendix  was  performed  by  Dr.  M.F.E.  Dillenius 
of  Nielsen  Engineering  & Research,  Inc. 
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Certain  output  has  been  added  to  allow  calculation  of  the  distri- 
bution of  lift  carryover  onto  the  body  and,  specifically,  to  aid  in  the 
determination  of  the  Kg  and  Xg  factors  described  in  the  main  text 
of  this  report.  Loading  information  for  each  body  ring  has  been  added 
to  the  section  of  the  code  concerned  with  the  properties  of  the  body 
influence  panels.  This  information  includes;  (1)  the  contribution 
from  each  ring  to  the  force  in  the  z-direction  (i.e.,  along  fin  1),  CZ, 
(2)  the  force  in  the  y-direction  (i.e.,  along  fin  4),  cy,  and  (3)  the 
moments  about  the  z and  y axes  taken  with  the  root  chord  leading  edge 
as  moment  center.*  The  moments  are  designated  CZMOM  and  CYMOM 
respectively.  In  addition,  the  accumulated  normal  and  side  forces  due 
to  lift  carryover  (CNB  and  CYB)  and  the  locations  of  the  axial  centers 
of  pressure  (XCP,CZ  and  XCP,CY)  are  given  as  a function  of  the  number 
of  rings  in  the  axial  direction.  To  determine  lift  carryover,  the  user 
should  look  at  the  contribution  to  CZ  and  CY  from  each  ring.  In  general, 
these  quantities  will  increase  for  each  successive  ring  in  the  axial 
direction  to  a maximum  and  then  decrease.  The  ring  experiencing  maximum 
CZ  and/or  CY  is  then  the  last  ring  over  which  lift  carryover  is  dominant 
(in  successive  rings  wake  effects  become  dominant) , The  lift  carryover 
is  then  determined  from  the  accumulated  CZ  and  CY  ever  the  affected  body 
rings  by  vectorial  addition  and  XCP,CZ  and  XCP,CY  are  the  corresponding 
centers  of  pressure. 

/ Z,CZ,CNB 


X y,CY,CYB 


Sketch  F.l.  Coordinate  system  used 
for  body  loading  due  to  lift. 


*Note  that  the  x,y,z  coordinate  system  used  here  is  for  rolled  body 
coordinates.  See  sketch  F.l. 
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F.3.  LOADING  CONVERGENCE  STUDY 
F.3.1  Modeling  of  Missile  Body 

A study  was  undertaken  to  investigate  the  effect  on  the  fin  loading 
of  the  number  of  line  sources/sinks  and  doublets  distributed  along  the 
body  centerline  to  model  the  missile  body.  Specifically,  the  magnitude 
of  the  velocity  component  normal  to  the  fin  surfaces  was  calculated  for 
different  numbers  of  body  singularities.  It  was  found  that  for  cross- 
flow  planes  aft  of  the  shoulder  of  an  ogive -cylinder  at  angle  of  attack, 
the  computed  upwash  varied  smoothly  with  distance  away  from  the  body. 

At  every  point,  the  magnitude  of  the  upwash  appeared  to  converge  rapidly 
with  increasing  numbers  of  body  singularities.  As  a rule  of  thumb,  the 
number  of  body  singularities  over  the  length  spanned  by  the  fins  should 
be  at  least  equal  to  the  number  of  constant  u-velocity  panels  laid  in  a 
chordwise  row  on  the  fins. 

F.3.2  Modeling  of  Fins  for  Wing-Alone  Case 

The  effect  of  the  number  of  constant  u-velocity  panels  distributed 
over  the  fins  and  body  interference  shell  were  studied  in  a systematic 
manner.  First  a wing-alone  was  considered.  The  normal  force  was  computed 
for  many  combinations  of  chordwise  and  spanwise  numbers  of  panels.  It 
was  concluded  that  the  normal  force  approaches  a definite  limit  as  the 
nvunber  of  chordwise  rows  is  increased  while  the  nurrtoer  of  panels  in  each 
row  is  kept  constant.  For  typical  missile  fin  planforms,  a layout  of 
9 panels  spanwise  and  4 panels  chordwise  should  give  good  results.  This 
recommendation  is  based  both  on  the  findings  from  the  convergence  studies 
for  a wing  alone  and  on  comparison  with  experimental  data. 

F.3.3  Modeling  of  Fins  and  Body  interference 
Shell  for  Wing-Body  Case 

Next,  the  effect  of  the  number  of  panels  on  the  fins  and  interference 
shell  for  a cruciform  fin-body  combination  was  investigated.  Based  on 
the  results  of  meuiy  calculations  of  the  normal  forces  acting  on  the  fins, 
it  was  concluded  that  the  size  of  the  panels  in  the  interference  shell 
should  be  comparable  to  the  size  of  the  panels  on  the  fins  in  the  chord- 
wise  row  at  the  fin-body  junction.  Thus,  the  number  of  panels  in  the 


interference  shell  in  the  eixial  direction  should  equal  the  number  of 
panels  in  a chordwise  row  on  the  fins.  The  number  of  panels  on  the 
circumference  of  the  interference  shell  should  be  selected  such  that 
the  panel  width  is  comparable  to  the  span  of  the  chordwise  row  of 
panels  on  the  fin  at  the  fin-body  junction. 

For  the  case  of  an  unrolled,  cruciform  fin-body  configuration  at 
angle  of  attack  and  with  the  vertical  panels  deflected  for  yaw  control, 
the  contribution  to  rolling  moment  from  the  deflected  fins  (termed 
direct  roll)  can  be  determined  with  good  accuracy  using  12  panels  on  the 
circvunference  and  4 along  the  body  length  spanned  by  the  fins,  in 
accordance  with  the  above  recommendations,  each  fin  should  be  covered 
with  3 or  4 chordwise  rows  with  4 panels  in  each  row.  In  contrast,  the 
contribution  to  the  rolling  moment  from  the  undeflected,  horizontal  fins 
(termed  reverse  roll)  was  very  small  for  the  case  considered.  Furthermore, 
the  convergence  was  not  monotonic.  In  this  instance,  the  number  of  panels 
required  to  reach  a limiting  value  exceeded  the  maximum  number  allowed 
by  the  computer  progrcun. 

F.4.  PRCX3RAM  LISTING 


A listing  is  given  below  of  the  two  subroutines  affected  by  the 
modifications  discussed  above.  Other  subroutines  were  changed  only  to 
make  COMMON  statements  compatible  or  to  remove  calls  to  the  system- 
supplied  subroutine  REQFL  and  to  place  the  aerodynamic  coefficient 
matrix  FVN  in  large  core  memory.  The  following  subroutines  have  b->en 
written  in  the  FORTRAN  IV  language  (029  punch)  for  the  CDC  7600  computer. 
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LIST  OF  SYMBOLS 


left  afterbody  vortex 

right  afterbody  vortex 

IR 

wing-alone  aspect  ratio 

\ 

see  equation  (B.29) 

a 

body  radius 

a 

e 

proportionality  coefficient  for  elliptic  span  loading 
equation  (C.2) 

; see 

proportionality  coefficient  for  linear  span  loading; 
equation  (C.l) 

see 

eucial-force  coefficient;  axial  force/q^^^Sj^ 

drag  coefficient;  drag/q^Sj^ 

^L 

lift  coefficient;  lift/q^Sj^ 

rolling-moment  coefficient;  rolling  moment/ q^^Sj^i^ 

pitching-moment  coefficient;  pitching  moment/q^^Sj^ij. 

S 

normal-force  coefficient;  normal  force/q^^Sj^ 

yawing-moment  coefficient;  yawing  moment/c^Sj^ij. 

(W) 

normal-force  coefficient  for  body  in  the  presence  of 
normal  force/q^^Sj^ 

fins; 

normal-force  coefficient  of  "wing  alone"  formed  by  joining 
two  opposing  fins;  normal  force/ q^^^Sj^ 

®*W(B) 

normal-force  coefficient  for  two  opposing  fins  in  the 
of  a body;  normal  force/<^Sj^ 

presence 

"«a 

slope  of  nose  noirmal-force  coefficient  at  = 0;  reference 

area  is  the  base  area  of  the  nose 

side-force  coefficient;  side  force/q^S^^ 

^6 

moderate-aspect-ratio  canard  fins  used  in  MICOM-MWC  tests; 
see  reference  4 and  figure  3(a) 

C7 

low-aspect-ratio  canard  fins  used  in  MIC(^-NWC  tests; 
reference  4 and  figure  3(b) 

see 
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LIST  OF  SYMBOLS  (CONTINUED) 


CNCl,CNC2,etc. 

CNi 

CNTl,CNT2,etc. 

CRMi 

CRMCl,CRMC2,etc. 

CRMTl,CRMT2,etc. 

Cl,C2,etc. 


normal-force  coefficients  associated  with  fins  cl, 

C2,  etc. 

normal-force  coefficient  for  fin  i 

normal-force  coefficients  associated  with  fins  Tl, 

T2,  etc. 

rolling -moment  coefficient  for  fin  i 

rolling-moment  coefficients  associated  with  fins  Cl, 
C2,  etc. 

rolling-moment  coefficients  associated  with  fins  Tl, 
T2,  etc. 

canard  fins;  also  vortices  associated  with  canard  fins 
two-dimensional  crossflow  drag  coefficient 


c 

r 

"=<>L 

“V 

K* 


K 


Op 


length  of  fin  root  chord 

elliptic  span  loading;  see  equation  (C.2) 

linear  span  loading;  see  equation  (C.l) 

span  loading;  see  equation  (B.5) 

span  loading;  see  equation  (B.14) 

crossflow  drag  force  on  an  elemental  length  of 
cylinder;  see  equation  (D.5) 

fraction  of  fin  leading-edge  or  side-edge  suction  which 
is  converted  to  vortex  lift 

wing-body  interference  factor  for  body  normal  force  for 

6 = 0,  / 0 

wing-body  interference  factor  for  fin  normal  force  for 

6 = 0,  0 

vortex  lift  parameter;  see  equation  (34) 

wing-body  interference  factor  for  bank,  5 = 0,  a ^ n 

c ' 

wing-body  interference  factor  for  body  normal  force 
for  6 / 0,  aj,  = 0 

wing-body  interference  factor  for  fin  normal  force 
for  6 ^ 0,  Oc  = 0 

see  equation  (B.27) 
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LIST  OF  SYMBOLS  (CONTINUED) 


t.e. 

MS  0 
N 

NV(x) 


‘Joo 

R 

r(x) 

Sc 

Si 

Sr 


leading  edge 
reference  length 
crossflow  Mach  mimber,  M sin  a 
rolling  moment 

free-strecun  Mach  number 

missile  station  zero;  origin  of  coordinate  system 
normal  force 

normal  force  due  to  viscous  crossflow 

number  of  vortices  in  the  crossflow  plane  located  at  cixial 
position  X 

difference  between  pressure  coefficient  on  suction  side 
of  fin  (or  body)  and  pressure  coefficient  on  impact  side 
of  fin  (or  body) 

dynamic  pressure  of  free  stream 

see  equation  (B.24) 

radius  of  circle  in  transform  plane;  see  equation  (A. 2) 

Vy^  + z^ 

radius  of  nose  at  axial  position  x 
body  planform  area  in  finned  section 

crossflow  area  upstream  of  first  set  of  fins 
pleuiform  area  of  fin  i 
reference  area 


Sjij  semispan  of  fin  on  body 

s(x)  fin  semispan  at  eucial  location  x 

Tl,T2,etc.  tail  fins 

T low-aspect -ratio  tail  fins  used  in  MICOM-NWC  tests;  see 

^ reference  4 and  figure  3(d) 


T 


14 


low-aspect-ratio  tail  fins  used  in  MICOM  tests;  see  refer- 
ence 16  emd  figure  3(c) 


w 


‘is 


t 

t.e. 


max 


V, 


N 


V 

w 

x,y,z 


XCTE 

XTLE 

^SlL 


MC 


linear 
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LIST  OF  SYMBOLS  (CCWTINUED) 


low-aspect-ratio  tail  fins  used  in  MICOM  tests;  see 
reference  16  and  figure  3(e) 


fin  spanwise  coordinate 
trailing  edge 

maximum  normal  velocity  which  can  be  induced  on  a fin 
volume  of  nose  from  tip  to  shoulder 
free-stream  velocity 


velocity  component  in  y-direction 
velocity  component  in  z-direction 


missile  body  axes;  x measured  positive  downstream  along 
body  rotational  axis,  y measured  positive  to  right  in  the 
plane  of  C4,  and  z measured  positive  upward  in  the  plane 
of  Cl;  origin  is  located  at  MS  0 


special  set  of  x,y,z  axes  for  0=0;  also  called  unrolled 
body  coordinates 


axial  position  of  the  canard  trailing  edge 

axial  position  of  the  leading  edge  of  the  tail  root  chord 

cucial  position  of  fin  hinge  line  measured  from  MS  0 


cixial  location  of  reference  for  moments  acting  on  missile 
measured  from  MS  0 


length  of  nose  from  tip  to  shoulder 

axial  location  of  initial  separation  of  body  nose  vorticity 
axial  position  of  center  of  pressure 


axial  position  of  center  of  pressure  for  body  normal  force 
due  to  presence  of  fins  measured  from  MS  0 


axial  position  of  the  center  of  pressure  for  fin  i measured 
from  MS  0 


axial  position  of  center  of  pressure  for  nose  normal  force 
due  to  potential  flow 


axial  location  of  center  of  pressure  for  nose  normal  force 
due  to  viscous  flow 


side  force 
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LIST  OF  SYMBOLS  (CC»}TINUED) 


side  force  due  to  viscous  crossflow 
coordinates  of  vortex  in  the  crossflow  plane 


spcuiwise  location  of  centroid  of  fin  trailing  vorticity 


y . or  z. 


lateral  location  of  the  center  of  pressure  for  fin  i measured 
from  body  centerline 


moment  arm  for  normal  force  induced  by  vortices 


fin  cuigle  of  attack 


included  angle  of  attack;  angle  between  x cucis  and  free- 
streeun  velocity 


”cr 


resulteuit  angle  of  crossflow  at  body  centerline  due  to  free- 
strecun  flow  and  external  vortices  exclusion  of  afterbody 
vortices 


ot- 


eq 


equivalent  angle  of  attack;  that  angle  of  attack  of  the  wing 
alone  for  which  its  normal  force  is  twice  that  of  the  fin; 
the  wing  alone  is  formed  by  joining  two  opposing  fins 
together 


a 


eq,P 

av^i(x.t) 

e 

Y 

r 

r’ 


equivalent  angle  of  attack  when  no  vortices  are  present 


local  angle  of  attack  induced  at  x,t  on  fin  i by  presence 
of  vortices 


average  of  a . (x,t)  over  chord 

fin  angle  of  sideslip 
see  equation  (B.26) 
vortex  strength 

nondimensional  vortex  strength;  vortex  strength/2T7V  a 


chamge  in  a^q  due  to  a collection  of  external  vortices  and 
their  images^ 


6 .6  .etc. 
1 ’ 2 * 


incidence  angles  associated  with  fins  Cl,  C2,  etc.,  for 
0 « 0,  6i  and  63  are  positive  for  trailing  edges  to  the 
right  viewed  from  the  rear,  and  63  and  6^  are  positive 
for  trailing  edges  down 


n 


nondimensional  spanwise  coordinate;  see  equation  (C.l) 


nondimensional  spanwise  position  of  fin  center  of  pressure; 
see  equation  (C.IO) 


235 


1 


w 


LIST  OF  SYMBOLS  (CONCLUDED) 


9 


A 


X 


o 


M. 


V 

Poo 

J 

J ■ 


<t> 


(JU 


direction  of  resultant  crossflow  velocity  at  body  centerline 
due  to  free-stream  flow  and  external  vortices  exclusive  of 
afterbody  vortices;  measured  counterclockwise  from  position 
y^-2ucis;  also  see  equation  (B.25) 

fin  leading-edge  sweep 

ratio  of  fin  tip  chord  to  root  chord;  also  complex  position 
in  the  transform  plane 

interactioii  coefficient  for  fin  rolling  moment  due  to  presence 
of  vortices;  see  equation  (B.18) 

complex  position  in  the  transform  plane 

density  of  fluid  in  free  stream 

complex  position  in  the  crossflow  plane 

complex  distance  between  a vortex  and  its  image  inside  a 
cylinder  nondimensionalized  by  the  cylinder  radius 

roll  angle,  angle  between  z axis  and  Zq  axis;  positive 
measured  clockwise  viewed  from  rear 

rnterdigitation  angle,  angle  between  canard  fin  Cl  and  tail 
fin  Tl  measured  clockwise  from  fin  Cl 

angular  rate  of  rotation  of  wing-body  combination 


Subscripts 

body  contribution  to  total  forces  and  moments  due  to  loading  seen 

by  body  in  finned  section 

d direct  flow 

fins  contribution  to  total  forces  and  moments  due  to  loading  seen 

by  all  four  fins  in  a finned  section 

L left  afterbody  vortex;  in  Appendix  C,  refers  to  linear  span 

loading 

R right  afterbody  vortex 

r reverse  flow 
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